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ABSTRACT 


This thesis deals with two different software packages to obtain the aerodynamic 
and control derivatives for a generic unmanned air vehicle (UAV). These data has a dual 
application. Firstly, it is required in the Mathworks’ Simulink 6-degree-of-freedom model 
of a generic unmanned air vehicle to develop a robust controller and do a variety of trade- 
offs. Secondly, is also needed to tune the parameters of the existing real-time controllers 
such as a Piccolo autopilot. 

The first approach explored in this thesis involves using the LinAir software 
program developed about a decade ago at Stanford University, the second one relies on 
the Athena Vortex Lattice package developed at Massachusetts Institute of Technology. 
The thesis applies two aforementioned packages to generate the aerodynamic data for two 
different-size UAVs, SIG Rascal and Thorpe Seeop P10B, emphasizing advantages and 
pitfalls of each approach, and further compares the obtained data with that of some other 
UAVs such as BAI Aerosystems Tern and Advanced Ceramics Corp. Silver Fox. The 


thesis ends with some computer simulations based on the obtained aerodynamic data. 
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I. INTRODUCTION 


A. BACKGROUND 


UAV modeling has been part of the design and modification process for new 
UAVs and has increasingly been used for rapid testing and verification. Changes and 
modifications are made to the model, evaluating it for intended function against the 
requirement. Modeling allows a faster verification and iteration of the design change 
cycle without having to conduct a flight test with a prototype UAV. It also reduces the 


cost incurred for each design change. 


At the Naval Postgraduate School’s Center for Autonomous Unmanned Vehicle 
Research, the two types of UAV simulations most widely used are the Matlab Simulink 


UAV model and the Piccolo Autopilot simulation. 


The Matlab UAV Simulink model is a generic 6 degree of freedom (6DOF) 
model for UAVs. This model is adopted in many research areas for UAV modification or 
operation scenario simulations. Figure 1. shows the top level model of the generic UAV 


model. 
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Figure 1. In-house UAV Matlab Simulink Model 


The test-beds used most often by the Center for Autonomous Unmanned Vehicle 
Research is the Rascal UAV, which uses the Piccolo Autopilot module. This autopilot 
module comes with a Piccolo simulator that allows the software to perform offline 


simulations. 





Figure 2. Piccolo IT (CCT Part No. 900-90010-00) 


UAV aerodynamic derivatives are required for both the UAV Matlab Simulink 


model and the Piccolo Autopilot simulator. UAV derivatives constitute a large portion of 


the input file for the Matlab Simulink 6DOF model, as shown in Appendix A. The same 


derivatives are required for the Piccolo Autopilot simulator. 


Various research projects have been carried out to determine the aerodynamic 


properties of the UAV. Figure 1. shows an extract of some UAV aerodynamic properties 


[1]. 











Pioneer” Bluebird*” FROG” Old Silver Fox gee 
Analytical approach Panel code 

Cio CLO _ lift coefficient at a= 0 0.385 0.02345 0.4295 0.3260 0.3280 0.3243 
C,, CLalpha lift curve slope 4.78 4.1417 4.3034 4.6800 5.0970 6.0204 
C,, CLa_dot lift due to angle of attack rate 2.42 1.5787 1.3877 0.8610 1.9300 1.9300 
Cz CLq _ lift due to pitch rate 8.05 3.9173 3.3) 2.5300 6.0300 6.0713 
Crs, CLDe _ lift due to elevator 0.401 0.4130 0.3914 0.3510 0.7380 0.9128 
Copp CDO = drag coefficient at C, = 0 0.06 0.0311 0.0499 0.0187 0.0191 0.0251 
A, Al drag curve coefficient at C; 0.43 0.1370 0.23 0.0000 0.0000 -0.0241 
A, A2 drag curve coefficient at C7 0.0413 0.0377 0.0692 
Cos, CDDe drag due to elevator 0.018 0.0650 0.0.0676 0.0486 0.1040 0.1000 
Cys CYb _ side force due to sideslip -0.819 -0.3100 -0.3100 -0.3100 -0.2040 = -0.3928 
Cys CYDr _ side force due to rudder 0.191 0.0697 0.0926 0.0613 0.1120 0.1982 
Cis, Clb dihedral effect -0.023 -0.0330 -0.0509 -0.0173 -0.0598 -0.0113 
Cr Clp roll damping -0.45 -0.3579 -0.3702 -0.3630 -0.3630 = -1.2217 
Cy Clr roll due to yaw rate 0.265 0.0755 0.1119 0.0839 0.0886 0.0150 
Cie, ClDa __ roll control power 0.161 0.2652 0.1810 0.2650 0.2650 0.3436 
Cis, ClDr — roll due to rudder -0.00229 0.0028 0.0036 0.0027 0.0064 0.0076 
Co  Cm0 pitch moment ata =0 0.194 0.0364 0.051 0.0438 0.0080 0.0272 
C,, Cma a due toangle of 519 4.0636 --0.5565—-0.8360---2.0510-1.9554 
Crag Cma_dot ? si —. seit ade S| -4.6790 -3.7115  -2.0900 —-5.2860_~—--5.2860 
C,; Cmq lie moment due to pitch 566 11.6918 -8.8818 6.1300 -16.5200  -9.5805 
C,,3, CmDe _ pitch control power 1.76 -1.2242 -1.0469 0.8490 ~—--2.0210-~—-2.4808 
Cup Cnb _ weathercock stability 0.109 0.0484 0.0575 0.0278 0.0562 0.0804 
Co Cnp = adverse yaw -0.11 -0.0358 -0.0537 -0.0407 -0.0407 -0.0557 
C,- Cnr yaw damping -0.2 -0.0526 = -0.0669 -0.0232 -0.0439 = -0.1422 
Cre, CnDa — aileron adverse yaw -0.02 -0.0258 -0.0272 -0.0294 -0.0296 -0.0165 
C,,, CnDr yaw control power -0.0917 -0.0326 -0.0388 -0.0186 -0.0377 -0.0598 





Table 1. | Non-dimensional aerodynamic and control derivatives for several UAVs (From 


[1].) 


B. OBJECTIVES 


The objective of this thesis is to determine the aerodynamics derivative using two 


software tools, LinAir and Athena Vortex Lattice (AVL). 
1. LinAir 


LinAir is a program that computes aerodynamic properties based on the model 
created. It is capable of generating the effect of angle of attack (AOA) and side slip angle 


on the force, drag coefficient, lift coefficient, etc. 


LinAir was developed first in 1982 and has been modified over the following 
years for various applications. It has been run on computers ranging from small laptops to 
VAXes and Crays. LinAir is now used in courses at several universities, by companies 
such as Boeing, AeroVironment, Northrop, and Lockheed, and by researchers at NASA 


to obtain a quick look at new design concepts [2]. 
2. Athena Vortex Lattice (AVL) 


Athena Vortex Lattice (AVL) is the other software program that was evaluated to 
compute the Rascal UAV aerodynamics derivative. AVL is a freeware created by Prof. 
Mark Drela of MIT. However, AVL alone is similar to LinAir, using a test file as the 
input source. It is possible to format and export the Rascal model in MS Excel for input 


into the AVL. 


However, the program AVL Editor, created by CloudCap Technology, serves as 
the graphic user interface to provide a window-based editor with manual coordinate input 
windows. The AVL Editor is able to provide visual inspection of the model created. The 


AVL software together with the AVL Editor can be _ found at 
http://www.cloudcaptech.com/resources_autopilots.shtm#downloads. 


C. REPORT STRUCTURE 


Section II discusses the modeling process of the P-10B model including LinAir 


and AVL. Section III presents the Rascal model including LinAir and AVL. Section IV 
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presents the results of simulink model responses to the data of LinAir and AVL. Section 
V reviews the limitation and constraints of LinAir and AVL, follow by conclusions. 


Finally, Section VI recommends directions for possible future research. 
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A. 


Il. MODELING OF P-10B UAV 


BACKGROUND 


The Thorpe Seeop P-10B is a Remotely Piloted Testbed (RPT), pictured in Figure 


3. An upcoming research program is planned for this UAV, and this thesis, which 


determines the aerodynamics derivative, has laid the groundwork for the future research 


work. 





Figure 3. Thorpe Seeop P-10B (From [3]) 


The base model from Thorpe Seeop Corporation is configurable to meet various 


needs of operation. Following is the specification of a configuration that will be modeled. 


[3] 


Model: Thorpe Seeop P-10B (http://www.seeop.com/) 
Wingspan: 21.25 feet 

Engine: 22 HP, two-cylinder, two-stroke 

Fuel cap: 3 US gal 

Max endurance: 2 hrs 

Speed: 70 mph 

Empty wt: 160 lbs 


Max takeoff w: 250 lbs 
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o Glimits: +3, -1.0 g’s 
o Payload weight: 50 to 85 Ibs 
o Payload volume: 11 x 11 x 23 inches 


o Runway requirement: 400 ft minimum for landing, 300 ft for 


takeoff 
B. GEOMETRY 


The graphical model of this UAV is provided in Appendix B. The airframe has a 
center of gravity at 7.5” from the leading edge of the wings, along the X axis, and from 


the center of the fuselage height. 


This UAV uses a Clark Y airfoil for the wing and a NACA 0008 for its horizontal 
stabilizer and vertical stabilizer. The Clark Y airfoil is a general airfoil that is widely 


used, and its. profile can be found” at  http://www.ae.uiuc.edu/m- 





selig/ads/coord_database.html. 

The profile was designed in 1922 by Virginius E. Clark. The airfoil has a 
thickness of 11.7 percent and is flat on the lower surface from 30 percent of chord back. 
The flat bottom simplifies angle measurements on propellers, and makes for easy 
construction of wings on a flat surface[4]. The airfoil profile data in Appendix E could be 


imported into the AVL easily for modeling. 


With the airfoil identified, the drag coefficient of the airfoil that consists of the 


following 2 portions can be identified. 


1. CDyin/ CDo - minimum drag that is achievable when the airfoil is in zero angle of 
attack (AOA). 
2. CDi-— induced drag as a result of increasing AOA and velocity of the airfoil. 


A typical drag coefficient model consists of CDO, CD1 and CD2. 
Cy =Cyy t+ AC, + A.C; 


These are also the input values for the LinAir model, matching the following modern 


drag equation from NASA [5] 


Cy = Cp min +(ax AR xE)Cr 
AR = WingSpan’ / WingArea 
é = efficiency = 0.9 


Therfore 
Coo = Cp min 


1 
mx ARxXé 


A, = 


254° | 
7620 — 
Assuming ¢=0.9 

: 1 

1x 8.467 x 0.9 





AR 8.467 


= 0.04177 


2 


From the above calculation, A2 = 0.04177. CDO is found to be about 0.01 based 


on Figure 4. 
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CLARK Y 
Figure 4. Clark Y Airfoil Chart (From [6]) 


C. MS EXCEL MODEL 


Data from the geometry of the UAV is used to create the MS Excel model. 
English units were adopted for the model, but it should be noted that the LinAir model is 
independent of units as long as the units are consistent within the model. A basic model 


was created with 22 elements that consist of the following: 
o 3elements for vertical fuselage 
o 3elements for horizontal fuselage 


o 2elements for left wing 
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o 2elements for right wing 

o 1 element for left aileron 

o 1element for right aileron 

o 2elements for horizontal stabilizer 
o 2 elements for left elevator 

o 2elements for right elevator 

o 2 elements for vertical stabilizer 

o 2 elements for rudder 


The reference area (Sref) is the projected wings area on the X-Y plane, is 
computed to be 7620 square inches. From the geometry model, the reference wing span 
(Bref) is determined to be 254 inches. Once the base model is completed, it is modified to 
adjust the coordination of the control surfaces (aileron, elevator, and rudder) with the 


input of angle of deflection. 
D. LINAIR MODEL 


Verification of the model is performed by exporting the model into a test file 
(space delimited format) and loading it with LinAir to check for visual error. Figure 5. 
shows the P-10B LinAir model; for the P-10B LinAir model input file, refer to Appendix 
C. 
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Figure 5. P-10B LinAir Model 


Once the model is completed, it is run for Alpha range from -10° to 15° in steps 
of 1° increments. Figure 6. shows the force contribution of each element and panel of the 
entire model. Most of the force contribution must be positive, especially for the wings, as 
it is the main lift contribution for the entire airframe. For a symmetrical airframe, the 
force contribution should be symmetrical as well. The only time it is asymmetric is when 
the aileron control surface is deflected, and the left side deflection is opposite to the right 


side. 
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Figure 6. Force contribution of entire P-10B model 


The next property to check is the relationship of Alpha (AOA) against lift 
coefficient (CL), as shown in Figure 7. The curve must have a positive gradient such that 


CL increases as AOA increases. 
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Configuration CL 
2 














Alpha (deg) 


Figure 7. P-10B LinAir Alpha vs. Lift coefficient curve 


The third check will be the drag coefficient (CD) vs. lift coefficient (CL). Figure 
8. shows a typical CD vs. CL polynomial curve. 
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Configuration CD 
0.15 





0.10 BP 


0.05 











0.00, 
Configuration CL 


Figure 8. P-10B LinAir Drag coefficient vs. Lift Coefficient curve 
Once all the checks were completed, the following runs were performed to gather 


the necessary data for coefficient analysis. 


1. Basic model with AOA varies from -10° to 15° and Beta varies from -10° 


to 15° at 1° step interval. 


2. Basic model with phat= 1 rad/s and Beta varies from -10° to 15° at 1° step 


interval. 


3. Basic model with qhat = 1 rad/s and AOA varies from -10° to 15° at 1° 


step interval. 


4. Basic model with rhat = | rad/s and Beta varies from -10° to 15° at 1° step 


interval. 


5. Deflection of Aileron for 5°, 10°, and 15° deflection and Beta varies from 


-10° to 15° at 1° step interval. 
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6. Deflection of Elevator for 5°, 10°, and 15° deflection and AOA varies 


from -10° to 15° at 1° step interval. 


7. Deflection of Rudder for 5°, 10°, and 15° deflection and Beta varies from - 


10° to 15° at 1° step interval. 


All data generated were imported into MS Excel for analysis, and the following 


coefficients were obtained. 


1. CLO is the lift coefficient, and at alpha is zero. This is the y-intercept of 
the Lift coefficient (CL) vs. AOA, as shown in Figure 9. 


2. CLalpha is the gradient of the Lift coefficient (CL) vs. AOA curve. This is 


a positive gradient curve, as shown in Figure 9. 















—e CL 
—Linear( CL ) 








CL 








10 15 20 





y = 0.088x - 0.0044 

















AOA 











Figure 9. Curve for CLO and CLalpha 


3. CLa_dot is the lift due to angle of attack rate. This value may be obtained 


by averaging the increase in lift due to AOA from -10° to 15°. 


4. CLq is the lift due to pitch rate. This is obtained by applying 1 rad/s pitch 
rate to the model without any control surfaces deflection. A curve of CL 
vs. AOA was plotted. The y-intercept is the lift due to pitch rate of the 
UAV, as shown in Figure 10. 
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y = 0.0857x + 0.0985 








5. CLDe is the lift due to elevator deflection. Various elevator deflection 
models were simulated, at 5°, 10°, and 15°. All three curves of CL vs. 
AOA were plotted in the same chart, and it is worth noting that they 
follow almost the same gradient. Differences of the y-intercept between 0° 
and 15° curve were obtained and divided by 15 (since total deflection is 


15°). This is the lift due to deflection of elevator per degree of deflection, 


Figure 10. 


as shown in Figure 11. 


Curve for CLq 
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Figure 11. 


Curve for CLDe 
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6. CDO0 is the drag coefficient at CL = 0. This coefficient is obtained from the 
CD vs. CL curve, which is typically a polynomial curve. CDO is the 
minimum drag coefficient when lift coefficient (CL) is at zero value, or y- 


intercept, as shown in Figure 12. 


7. Al is the drag curve coefficient at CL. This coefficient is obtained from 
the CD vs. CL curve. Al is the coefficient of first order lift coefficient, 


CD=A1*CL42+A2*CL+CDO, as shown in Figure 12. 


8. A2 is the drag curve coefficient at CL2. This coefficient is obtained from 
the CD vs. CL curve. A1 is the coefficient of second order lift coefficient, 


CD=A1*CL42+A2*CL+CDO, as shown in Figure 12. 


























CD 




















-1 -0.5 0 0.5 1 1.5 
CL sy = 0.0774x? + 0.0004x + 0.0129 











Figure 12. Curve for CDO, Al and A2 


9. CDDe is the drag due to elevator deflection. Various elevator deflection 
models were simulated, at 5°, 10°, and 15°. All three curves of CD vs. 
AOA were plotted in the same chart, andonce again they have almost the 
same gradient. Differences of the y-intercept between 0° and 15° curve 
were obtained and divided by 15 (since total deflection is 15°). This is the 
drag due to deflection of elevator per degree of deflection, as shown in 


Figure 13. 
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Figure 13. Curve for CDDe 


10. CYb is the side force due to sideslip. This coefficient can be obtained from 
the gradient of the CY vs. Beta curve. Typically the y-intercept is zero if 
the model is stable and no side force is expected for zero side slip, as 


shown in Figure 14. 
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Figure 14. Curve for CYb 


11. CYDr is the side force due to rudder deflection. Various rudder deflection 
models were simulated, at 5°, 10°, and 15°. All three curves of CY vs. 


Beta were plotted in the same chart and follow almost the same gradient. 
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Differences of the y-intercept between 0° and 15° curves were obtained 
and divided by 15 (since total deflection is 15°). The value obtained is the 


side force due to one degree of rudder deflection, as shown in Figure 15. 
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Figure 15. Curve for CYDr 


12. Clb is the dihedral effect. The dihedral effect coefficient could be obtained 
from the gradient of the CR vs. Beta curve. Typically this curve has zero 


y-intercept, as shown in Figure 16. 
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Figure 16. Curve for Clb 
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13. Clp is the roll damping. Roll damping is the roll effect of the UAV as a 
result of applying 1 rad/s roll rate without any control surfaces deflection. 
A polynomial curve of CR vs. Beta was plotted. The y-intercept is the Roll 
damping of the UAV, as shown in Figure 17. 





—e— CR 
—Poly.( CR ) 


CR 





Beta y = 1E-05x? - 0.002x - 0.0707 











Figure 17. Curve for Clp 


14. Clr is the roll due to yaw rate. This is obtained by applying 1 rad/s yaw 
rate and without any control surfaces deflection. A curve of CR vs. Beta 
was plotted. The y-intercept is the Roll damping of the UAV, as shown in 


Figure 18. 
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Figure 18. Curve for Clr 
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15. ClDa is the roll control power due to aileron deflection. Various aileron 
deflection models were simulated, at 5°, 10°, and 15°. All three curves of 
CR vs. Beta were plotted in the same chart and follow almost the same 
gradient. Differences of the y-intercept between 0° and 15° curves were 
obtained and divided by 15 (since total deflection is 15°). The value 
obtained is the side roll rate due to one degree of aileron deflection, as 


shown in Figure 19. 
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Figure 19. Curve for ClDa 


16. ClDr is the roll due to rudder deflection. Various rudder deflection models 
were simulated, at. 5°, 10°, and 15°. All three curves of CR vs. Beta were 
plotted in the same chart and follow almost the same gradient. Differences 
of the y-intercept between 0° and 15° curves were obtained and divided by 
15 (since total deflection is 15°). The value obtained is the side roll rate 


due to one degree of rudder deflection, as shown in Figure 20. 
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Figure 20. Curve for ClDr 


17. CMO is the pitch moment at AOA = 0. CMO is the y-intercept for the 
curve between moment coefficient (CM) and AOA. It is obtained from the 


run where there is no control surfaces deflection, as shown in Figure 21. 


18. CMa is the pitch moment due to angle of attack. It is the gradient of the 
curve between moment coefficient (CM) and AOA. Typically this is a 
negative gradient. It is obtained from the run where there is no control 


surfaces deflection, as shown in Figure 21. 
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Figure 21. Curve for CMO and CMa 
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19. CMa_dot is the pitch moment due to angle of attack rate. This value may 
be obtained by averaging the increase in pitching moment (CM) due to 


AOA from 0° to 15°. 


20. CMq is the pitch moment due to pitch rate. This is obtained by applying 1 
rad/s pitch rate to the model without any control surfaces deflection. A 
curve of CM vs. AOA was plotted. The y-intercept is the pitch moment 
due to pitch rate of the UAV, as shown in Figure 22. 
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Figure 22. Curve for CMq 


21. CMDe pitch control power is the moment due to elevator deflection. 
Various elevator deflection models were simulated, at 5°, 10°, and 15°. 
All three curves of CM vs. AOA were plotted in the same chart and follow 
almost the same gradient. Differences of the y-intercept between 0° and 
15° curve were obtained and divided by 15 (since total deflection is 15°). 


This is the lift due to deflection of elevator per degree of deflection, as 


shown in Figure 23. 
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Figure 23. Curve for CMDe 


22. CNb is the weathercock stability. It is the ability of an aircraft to return to 
its previous heading after being yawed as a result of wind effect [7]. This 
coefficient could be obtained from the gradient of the CN vs. Beta curve. 


This curve typically has a zero y-intercept, as shown in Figure 24. 
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Figure 24. Curve for CNb 


23. CNp is the adverse yaw of the UAV. This is obtained by applying | rad/s 
roll rate to the model without any control surfaces deflection. A curve of 
CN vs. Beta was plotted. The y-intercept is the adverse yaw effect of the 


UAV, as shown in Figure 25. 
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Figure 25. Curve for CNp 


24. CNr is the yaw damping of the UAV. This is obtained by applying | rad/s 
yaw rate of the model without any control surfaces deflection. A curve of 
CN vs. Beta was plotted. The y-intercept is the yaw damping of the UAV, 


as shown in Figure 26. 
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Figure 26. Curve for CNr 


25. CNDa is the aileron adverse yaw of the UAV. It is the aileron control 
power due to aileron deflection. Various aileron deflection models were 
simulated, at 5°, 10°, and 15°. All three curves of CN vs. Beta were 


plotted in the same chart and follow almost the same gradient. Differences 
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of the y-intercept between 0° and 15° curves were obtained and divided by 
15 (since total deflection is 15°). The value obtained is the yaw control 


power due to one degree of aileron deflection, as shown in Figure 27. 
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Figure 27. Curve for CNDa 


26. CNDr is the yaw control power of the UAV due to rudder deflection. 
Various rudder deflection models were simulated, at 5°, 10°, and 15°. All 
three curves of CN vs. Beta were plotted in the same chart and follow 
almost the same gradient. Differences of the y-intercept between 0° and 
15° curves were obtained and divided by 15 (since total deflection is 15°). 
The value obtained is the yaw control power due to one degree of rudder 


deflection, as shown in Figure 28. 
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Figure 28. Curve for CNDr 
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Pioneer”Bluebird*|FROG""| Old Silver Fox |N¢™ S#¥€"|_ pyop 
Nomenclatures Fox 


Cana, || PEE omni ne ee) 5 19.'|' 1 6636 | 0.9565 | 0836 | D051 | 19554 
attack ‘ -2.4694 

Cada |e eee sii -4.679 |-3.7115]} -2.09 | -5.286| -5.286 
attack rate -2.0695 

Cig | PP omen ake epee aes | 11.6098 |-8.8818| 6.13 | 1650 | 08808 
rate i -8.4569 


Table 2. Consolidation of P-10B coefficient Using LinAir (After [1]) 





E. AVL MODEL 


The first step in AVL modeling is to define the fuselage of the UAV. As AVL is 
only able to model circular fuselages, equivalent areas were used. In addition, the 
dimensions for the P-10B were scaled down by 50% since an AVL fuselage is limited to 
100 units or less. The following parameters were used to create the fuselage, as shown in 
Figure 29. 
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4 Create New Body 





Body Name | Rascal 
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Overall Length [D] 





Max Cross Section 4rea [B] 
Distance From LE [A] 
TE Surface Area [C] 
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San 





Figure 29. AVL GUI to create fuselage 


The next step was to define the airfoil required for the model. The airfoil editor 
was activated, as shown in Figure 30. Defining of NACA airfoil is done by simply 
entering the airfoil number. NACA 0008 is defined for vertical and horizontal stabilizers. 
The wing uses a Clark Y airfoil. Non-NACA airfoils can be defined by loading the airfoil 
profile that comes in the form of a text file. It is important to note that the airfoil profile 
should start from x = 1.00000 (trailing edge) and proceed towards zero value going 
around the leading edge and back to the trailing edge. The connecting point should be at 
the trailing edge, shown as a small circle in Figure 31. 


 Kirfoil Editor 


H/C yic Airfoil Flat Panel 
1 0.000000 0.000000 Preview 


ya) 1.000000 0.000000 








wh 





Figure 30. AVL Airfoil Editor 
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0.945628 0.011660 
0.937646 0.013110 
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Figure 31. AVL NACA Clark Y Airfoil Profile Input 





Once the airfoils are defined, the next step is to define the surfaces. This includes 
all surfaces, excluding the fuselage that was defined previously. The surfaces were 
defined using a surface editor, as shown in Figure 32. The coordinates from the MS Excel 


model were used as an input to the surface editor. 
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* Surface Editor 





Surface Rudder 


v ¥-Symmetric [_] 








Property 
=) Sections 
(=) Section 1 
- Position (Leading Edge} 
Chord 
Incidence Angle 
Airfoil 
&- Control Surface 1 
Name 
- Hinge Vector 
Chord Fraction 
Deflection 
Gain 
H- Section 2 
(4) Section 3 
4 Section 4 
=) Geometric Transformation 
fH Scale 
4) Translation 
Incidence Angle 
(=) Vortex Lattice 
Chord Yortices 
Chord Spacing 
Span Vortices 
Span Spacing 


ee ee er ee 





“an 


(0.00000, 0,00000, -3.00000) 
1.50000 

0.00000 

Flat Panel 


Rudder 

(0.00000, 0.00000, 0.00000) 
0.00100 

Symmetric 

0.00000 


(1.00000, 1.00000, 1.00000) 
(46.60000, 0.00000, 1.50000) 
0.00000 


12 
2.00000 
64 


0,00000 rc: 








Apply 


Delete Surface} | Delete Section | | Delete Control 





Figure 32. AVL GUI for surface editor 


The complete P-10B model in AVL is shown in Figure 33. 
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Figure 33. P-10B Model in AVL 


Unfortunately, AVL is unable to generate positive results for the P-10B model. 
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Hl. MODELING OF RASCAL UAV 


A. BACKGROUND 


Rascal UAV is an in-house integrated UAV that is often used as a test bed for 
various UAV operation concepts and scenarios. The Rascal UAV uses a remote 
controlled aircraft (Sig Rascal) as the airframe platform, and integrates the following 


components: 
o Cloud Cap Technology Piccolo II Autopilot Avionics 
o GPS 
© pitot-static probe 
O magnetometer 
o 900 MHz radio modems for GCS link 
o PC/104 onboard computers (2) 
o 2.4 GHz “Mesh” wireless networking card 
o Video camera with custom 2-axis gimbals 
o Pelco network digital video server 


The integrated Rascal UAV weighs 13 pounds, measures 6 feet long, has a 9 foot 
wingspan, and has an endurance of 2 hours. The Rascal UAV has been a test bed for 


various UAV capabilities research and testing. [8] 
B. GEOMETRY 


Since the Rascal UAV is a hobby remote controlled aircraft, and no engineering 
drawing can be found, physical measurement of the airframe is performed. The UAV is 
divided into fuselage, wings, horizontal stabilizer, and vertical stabilizer. The 
measurements are converted into the graphical model shown in Appendix B, using MS 
Visio. 
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It was also found that the center of gravity of the airframe is 2.5 inches from the 


leading edge of the wing, along the X axis and center of the airframe height. 


The wing uses NACA 2312 airfoil with a 2% camber at a distance of 0.3% of the 
total cord length from the leading edge, with a 12% thickness. The cross sectional view 


of the airfoil is as shown in Figure 34. 


Index Name > 1232 -A00 
Serial No. : NO1-0076-0 
Airfoil Name : N.A.C.4. 2312 
Airfoil Info. : Computed 


Figure 34. Cross sectional view of NACA 2312 airfoil[9] 


Similarly to the P-10B model, the next step is to determine the drag coefficient 
after identifying the airfoil number (reference Figure 35. CDO = 0.017 for 1° AOA since 


the wingspan is mounted at 1 degree angle. 
a 1 
xX ARXE 
_ 110.4 
“T1227 
Assuming ¢=0.9 
1 


= ——__————_ = ().020668746 
mw xX17.11171x0.9 


2 


=17.11171 
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Figure 35. NACA 2312 Drag Coefficient vs. Mach [10] 


Alternatively the CD parameters could be computed using the AVL program. A 
portion of the NACA airfoil was created in AVL program as shown in Figure 36. The 
CD vs. CL curve obtained from the results is shown in Figure 37. From the curve, the CD 


coefficients are as follows, 
o CDO=0.000033 
o CD1=-0.0043 


o CD2=0.0002 
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+ Untitled Aircraft - AVL Model Editor [D:/CourseWork/MEOB20 Thesis/Rascal/AVL/AVL Editor. 1.0. 1/AVLEditor/NACA231 2.avl] 
File Model View Help 
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Scale: [Custom | XZPlane Y-ZPlane X-Plane 


Figure 36. NACA 2312 airfoil in AVL 





CD vs CL for NACA 2312 




















CL y = 0.0002x? - 0.0043x + 0.00033 











Figure 37. CD vs. CL for NACA 2312 using AVL 
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C. MS EXCEL MODEL 


Once again, English units were adopted. The MS Excel model created has 28 
elements that consist of the following: 
o 5 elements for vertical fuselage 
o 4elements for horizontal fuselage 
o 3 elements for left wing 
o 3 elements for right wing 
o 1 element for left aileron 
o 1 element for right aileron 
o 3 elements for horizontal stabilizer 
o 2 elements for left elevator 
o 2 elements for right elevator 
o 2 elements for vertical stabilizer 


o 2 elements for rudder 


The Sref (wings area) is computed to be 712.27 square inches. From the geometry 
model, the Bref (wing span) is determined to be 110.4 inches. Once the base model is 
completed, it is modified to adjust the coordination of the control surfaces (aileron, 


elevator, and rudder) with the input of angle of deflection. 


D. LINAIR MODEL 


Figure 38. shows the LinAir model for the Rascal UAV, with the input file 
attached in Appendix C. The following three charts were generated to verify that the 


model is correct. 


o Contribution of force component in the Y axis for the entire model; 


refer to Figure 39. 
o Alpha vs. lift coefficient (CL); refer to Figure 40. 


o Drag coefficient (CD) vs. lift coefficient (CL); refer to Figure 41. 
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Figure 38. Rascal LinAir model 
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Figure 39. Force contribution along Y axis for Rascal LinAir Model 
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Figure 40. Lift coefficient vs. AOA for Rascal LinAir model 
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Figure 41. Drag coefficient vs. lift coefficient for Rascal LinAir model 


Following the same procedure as P-10B model in Section II, the following 


coefficients for Rascal were determined: 
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Abbreviations 
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4.78 
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8.05 

0.401 


0.018 
b sidesli, -0.819 
0.191 
b -0.023 


drag due to elevator 
side force due to sideslip 
side force due to rudder 


dihedral effect 


aN 
(=) 


Bluebird” 


0.02345 


4.1417 
1.5787 
3.9173 
0.413 
0.0311 
0.137 


0.065 
-0.31 
0.0697 
-0.033 


FROG”! 


Analytical approach 


0.4295 
4.3034 
1.3877 
3.35 
0.3914 
0.0499 
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0.0.0676 
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0.0926 
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Fox 
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is 
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Analytical approach Panel code 
-0.45 -0.3579 | -0.3702 | -0.363 : -0.363 | -1.2217 | -3.2258 
0.265 0.0755 | 0.1119 | 0.0839 : 0.0886 0.015 0.4011 
0.161 0.2652 0.181 | 0.265 | 0.265 0.3436 0.4771 


ClDr roll due to rudder -0.00229 | 0.0028 | 0.0036 | 0.0027 ' 0.0064 | 0.0076 | 0.0057 
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Table 3. | Consolidation of Rascal coefficients using LinAir (After [1]) 


es) 


: AVL MODEL 


The same procedures were performed as in the case of P-10B model for AVL. 
The fuselage was defined in full scale, since the entire UAV is less than 100 inches in 
length. NACA 2312 was defined while a flat plate was used for vertical and horizontal 


stabilizers. The complete Rascal model is shown in Figure 42. 
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Figure 42. Rascal Model in AVL 


Using the same procedure as set out in P-10B, the following coefficients for 


Rascal were determined: 
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Table 4. Consolidation of Rascal coefficients using AVL (After [1]) 
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IV. SIMULINK MODEL 


A. MOMENT OF INERTIA 


With the above coefficient computed from the LinAir model or AVL model. The 
next parameter required by the simulink 6DOF model is the three moment of inertia for 
the x axis, y axis, and z axis. As the information of the moment of inertia for both P10B 


and rascal are not available, an estimation of the moment of inertia were done using the 


following formulas, 


2 
mass span 





=I XX ,ref 2 
MASS, SPAN, 

— mass lenght* 
yy — *YyY,ref 2 
mass, lenght,,, 


2 
mass span 





[T, =f ZZ,ref 2 
MASS jy SPAN, 





Bluebird”| FROG" | Old Silver Fox| N€¥ Silver | prog | Rascal 


Fox 
S 1.6260 0.6360 
b_ |wingspan, m 3.2250 2.4100 
c 0.500 0.264 
m 30.7 10.0 
F 19 147 
i 27.0 26.00 
Is 17.0 1.00 
114 087 


Izz |yaw inertia, kgm” 27.1 953: 1.40 
Table 5. Moment of inertia for P1OB(After [1]) 
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Using the code as shown in Appendix G, the trim condition for each UAV are 


determined as shown in Table 6. 













Rascal Rascal 


LinAir 





Velocity (m/s) 
OA Trim (degree) 


[Elevation trim, detrim 
(degree) 
Thrust Trim (N) 
















Table 6. | Aerodynamic model parameter at trim condition 


The trim parameter for the Rascal UAV coefficient generated by AVL, is logical. 
The elevator has to deflect by -34 degree, but the model has limited the maximum 


deflection to +15 degree. In addition, the thrust required is an order of magnitude larger 


then the LinAir’s result. 
B. P10B SIMULATION RESULTS 


With the aerodynamic coefficient, moment of inertia, and trim condition 
parameter, simulation were performed to verify that the aerodynamic coefficient 


determined are good enough to be used for simulation. 
1. Trim condition 


The first simulation was to determine whether the UAV is stable at trim condition. 
Figure 43. shows that the longitudinal channel of the UAV. It demonstrated an initial 
short period of disturbance that dies of in about 1 second. Following on is a long period 
of small oscillation. The magnitude of this oscillation is small, less than 10 meter for 


altitude. Figure 44. shows that the at trim condition, the horizontal channel of the UAV is 


constant at zero. 
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Figure 44. P10B lateral channel at trim condition 
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2. Change in Thrust Parameter 


Once the trim condition of the UAV was verified, the thrust of the model is varies 
to verify the behavior of the model. Figure 45. shows that with a 1ON increase in thrust, 
the UAV start to climb in altitude. This is due to more lift generated with additional 


speed generated from the additional thrust. The rest of the vertical channel does not 


change much. 
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Figure 45. P10B longitudinal channel with 22% (10N) throttle increase 
3. Change in Elevator Deflection 


Figure 46. shows that with -1 degree elevator deflection, the UAV generates more 


lift. Thus the UAV’s altitude increases. However, the deflection generates more drag that 


causes the speed to drop slightly. 
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Figure 46. P10B longitudinal channel with 1 degree elevator deflection 
4. Change in Aileron Deflection 


Figure 47. shows that as aileron deflect, the UAV start to roll. The rolling of the 


UAV also causes the UAV to experience a certain yaw effect. 
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Figure 47. . P1OB lateral channel with | degree aileron deflection 


5. Change in Rudder Deflection 


As the rudder deflect, the UAV yaw as shown in Figure 48. . At the same time, 


the UAV also experience roll effect. 
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C. RASCAL SIMULATION RESULTS 
The simulation result of the Rascal model is very similar to the P10B. 
1; Trim Condition 


Figure 49. shows that the longitudinal channel of the UAV at trim condition. It 
demonstrated an initial short period of disturbance that dies of in about 1 second. 
Following on is a long period of small oscillation. The magnitude of this oscillation is 
small, less than 10 meter for altitude. Figure 50. shows that the at trim condition, the 


horizontal channel of the UAV is constant at zero. 
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2. Change in Thrust Parameter 


A IN additional thrust was introduced to the model. The UAV start to climb in 
altitude as a result of additional lift generated from the higher velocity, as shown in 


Figure 51. . The rest of the vertical channel does not change much. 






























































10 
co 
ax 5 
fo} 
< 
0 
0 10 20 30 
Time (s) 
50 @ 200 
a 
ce o 
© 
& oF \ 2 0 
2 © 
a S 
-50 = -200 
0 10 20 30. & 10 20 30 
Time (s) Time (s) 
350 ~ 40 
wn 
: e 
g 300 g 30 ’ 
= om 
< 250 D 20 
10 20 30 0 10 20 30 
Time (s) Time (s) 
Figure 51. Rascal longitudinal channel with 21% (1N) throttle increase 
3. Change in Elevator Deflection 


Figure 52. shows that with -1 degree elevator deflection, the UAV generates more 


lift as the elevator deflect downwards. Therefore, the altitude increase and the speed drop 


due to increase drag. 
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Figure 52. Rascal longitudinal channel with -1 degree elevator deflection 
4. Change in Aileron Deflection 


Figure 53. shows that as aileron deflect, the UAV start to roll. The rolling of the 


UAV also causes the UAV to experience a certain yaw effect. 
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Figure 53. -Rascal lateral channel with 0.1 degree aileron deflection 


5, Change in Rudder Deflection 


As the rudder deflect, the UAV yaw as shown in Figure 54. . At the same time, 


the UAV also experience roll effect. 
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V. DISCUSSION AND CONCLUSION 


Various limitations were found during the modeling of the UAV in LinAir and 


AVL software. Following are listed the constraints and limitations of LinAir and AVL. 


A. LINAIR LIMITATIONS AND CONSTRAINTS 


The following were observed during the modeling in LinAir: 


1. 


The number of panels for each element must be input into the system 
manually, but the number of panels cannot be determined before 
modeling. A “wrong” combination of numbers of panels for the elements 


will result in drastic errors in the output. 


Knowledge of the control keys is necessary to control the software, 1.e., 


“L” for moving left, “P” for pitch, etc. 
The current version used does not support long filenames. 


The software is unable to determine the drag vs. lift relationship for 
known airfoils. Conversion of airfoils into drag coefficients is required, 


and forms part of the input file. 


. When exporting the MS Excel model as an input file for LinAir, it was 


noted that the cell spacing of the Ms Excel worksheet has to be wide 
enough such that numbers would not be seen as joined together with the 
number in its neighbor cell. This is because the visual width of the MS 
Excel’s cell is equivalent to the spacing for the space delimited text file. 
Therefore, a cell full with numbers will appear merged with its neighbor 


cell after exporting. 


The model created in LinAir is dimension independent. Either SI units or 
English units should be used throughout the input parameter. However, the 


same unit must be consistently used throughout the entire model. 


my 


2 


10. 


11. 


12. 


Root to Tip direction must be in the positive direction of Y and Z axes for 
the horizontal and vertical plane respectively. i.e., the LE root is smaller 
than the LE tip. This means that the positive Y axis is toward the right side 
of the UAV (viewing from the rear of the UAV). Having an opposite sign 
convention will result in a negative force element while the CL vs. AOA 


remains normal. 


The Nelem field is important, as LinAir determines the total number of 


elements from this field. 
All elements of the same part should have the same drag coefficients. 


It is better and easier to troubleshoot if fewer panels are defined for each 


element. 


As far as possible, elements should not be stacked against each other along 


the X axis under the same surface, i.e., wing, vertical stabilizer, etc. 


Every loading of the model will refresh the data file. Running of the 
simulation continuously without reloading will result in the output file 
been congested. It is recommended that the model be reloaded and data 


extracted at each interval before proceeding to the next configuration. 


B. AVL LIMITATIONS AND CONSTRAINTS 


The following were observed during the modeling in AVL: 


ly 


AVL uses +X axes from the leading edge to the trailing edge of fuselage, 


and +Y axes from the left wing to the right wing, and +Z axis upwards. 


AVLs have only circular fuselages. Non-circular fuselages are 
approximated using circular fuselage measurements. AVL’s fuselage GUI 
is only able to create a fuselage with up to a length of 100 units. In order 
to overcome this constraint, the physical value is scaled down to below 
100. Subsequently, the actual scale can be restored by indicating an 


equivalent upscale on the X, Y and Z scale value. 
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3. Similar to LinAir, the leading edge and trailing edge direction must be the 


same for joint parts such as left wing and right wing. 


4. No drag information of the airfoil is required, as the system is able to 
accept airfoil geometry to determine the drag properties. Defining on 
NACA airfoil is done by simply entering the NACA number. AVL also 
allows import of other airfoil geometry by either manual entering for 


import form text file. 
5. The AVL requires a much longer duration for each simulation run. 


6. The total dimension for any axis must be less than 100 units. AVL GUI 
will output the data into a text file (refer to Appendix F) before activating 
the AVL program to obtain the data from the text file created. Having any 
dimension more than 100 units will result in this value merging with its 


neighboring value, and cause errors in the AVL program. 


7. AVL is unable to create control surface with sudden increase in cord 


fraction, such as the case for the P1OB and the rudder for Rascal. 


8. AVL requires high processing power and/ or time. 


C. CONCLUSION 


Unfortunately AVL is unable to generate result for P1OB and the coefficient 
generated for Rascal does not produce meaningful trim condition. As such, comparison 


between LinAir and AVL is not possible. 


LinAir generated coefficient for both P1OB and Rascal has generated satisfactory 


trim condition parameter, as well as simulation result. 


However, LinAir is preferred over AVL based on the experience of using the 
software. This is due to the significant different in speed of the simulation and graphic 


user interface of the software. 
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VI. RECOMMENDATION 


The following research area are recommended for future work on aerodynamics 


derivative, 


A. AVL 


This research has not produce satisfactory result from AVL for comparison 
between LinAir and AVL. It is recommended to further evaluate the AVL software and it 


integration with Piccolo simulator. 


B. VERIFICATION WITH ACTUAL DATA 


In order to verify that the 6DOF model with the aerodynamic and control 
derivative is a good representative of the actual model, a comparison of the simulated 


path with an actual flight data may be performed. 


C. FUTURE WORK ON P10B AND RASCAL 


The result from this thesis on P1OB and Rascal could be used in the controller 
design for these two UAV using Simluink model. The results could also be used in the 


Piccolo simulator to assess the flight quality of the UAV before actual flight. 
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APPENDIX A. MATLAB SIMULINK INPUT FILE 


























A. P-10B SIMULINK MODEL INPUT FILE 

GL; 

T=0.05; 

r_lim=7; 

%& Initial Conditions in ENU (all vector data is represented as a 
column vectors) 

Pos_0 = [-1000; -500.0; 300]'; % Initial position vector (m) 
Vel_0 = [30.; 0; O]'; % Initial velocity vector (m/s) 
Euler_O = [0; 0; 0O]'*pi/180; % Initial Euler angles (rad) 
Omega_0 = [0; 0; O]'; 6 Initial Omega 

(rad/s) 

POR_O = [Os Os O] 's % Initial Omega (rad/s) 
Vb_0 = [30; 0; O]'; % Initial body-velocity vector 
(m/s) 


& Mass and Geometric 








Parameters recomputation 
































iS) = 4.916119; % surface area of wing (m2) 

span = 6.4516; %& wingspan (m) 

chord = S/span; % chord (m) 

mass = 63; % gross weight (kg) 

Ixx = 45; %$ main moment of inertia around axis Ox (kg*sq.m) 
Iyy = 26; % main moment of inertia around axis Oy (kg*sq.m) 
zz = 61; %$ main moment of inertia around axis Oz (kg*sq.m) 
& Aerodynamic Derivatives (all per radian) 

CLO = 0.2521% & lift coefficient at a = 0 = 0.0003; 

CLa = 5.0420; % ft curve slope 

CLa_dot = 4.8890; %& lift due to angle of attack rate 
CLgq = 5.6436; % lift due to pitch rate 

CLDe = 0.3621; % lift due to elevator 

CDO 0.0129; % drag coefficient at a = 0 

Apolar = 0.0774; %& drag curve slope (A2) 

Al = 0.0004; 

CYb = -0.2922; % side force due to sideslip 

CYyDr = 0.1587; % sideforce due to rudder 

Clb = -0.1146; % dihedral effect =-0.0132 

Clp = —-4.0508; % roll damping 

Clr = 0.3495; %& roll due to yaw rate 

ClDa = 0.4509; ap ae) control power 

C1lDr = 0.0052; % rol due to rudder 

Cm0 = 0.0974; %& pitch moment at a = 0 =>0.0652 

Cma = -2.4694; % pitch moment due to angle of attack 
Cma_dot = -2.0695; % pitch moment due to angle of attack rate 
Cmq = -8.4569; %& pitch moment due to pitch rate 

CmDe = -1.0273; % pitch control power 

Cnb = 0.0573; %& weathercock stability = 0.075 

Cnp = -0.2693; %& adverse yaw 

Cnr = -0.3782; %& yaw damping 

CnDa = -0.0045; %& aileron adverse yaw 
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CnDr = -0.0586; 


CLDE£ = 0; % 
CmD £ = 0; 

% Standard Atmosphere 
ISA_lapse = .0065; 
ISA_hmax = 2000; 
ISA_R = 287; 
ISA_g = 9.815; 
ISA_rho0 = 1.225; 
ISA_PO = 101325; 
ISA_TO = 289; 


o2..8. 02. C..6. 20.0.0. 9.0. 0.0.9 0.0 2.00. 


° 
O 


22 


2.e 


yaw con 


F2AVDD 


F2RV9 9 9 


trol power 


Lapse rate degC/m) 
Altitude limit m) 
Gas Constant degk*m*m/s/s) 





( 

( 

( 
Gravity (m/s/s) 

( 

( 

( 


Density at sea level kg/m/m/m) 
Sea-level Pressure N/m/m) 
Sea-level Temperature (degk) 











YOVIAAAAAAAAAIAA DD 


9.9.0.0. 9. C0. 0.0.9.0. 0.0.0 9: 0..0. 























(onker 

T=0.05; 

r_lim=7; 

& Initial Conditions in ENU (all vector data is represented as a 
column vectors) 

Pos_0 = [-1000; -500.0; 300]'; % Initial position vector (m) 
SVel_0 = [25.7 O¢ 0] "s & Initial velocity vector (m/s) 
EBuler_O = [0; 0; O]'*pi/180; % Initial Euler angles (rad) 
Omega_0 = [0; 0; O]'; 6 Initial Omega 

(rad/s) 

POQR_O = [Os Of O] 's % Initial Omega (rad/s) 
Vb_0O = [27; 0; O]'; % Initial body-velocity vector 
(m/s) 


& Mass and Geometric Parameters recomputation 


s = 0.459528; 
span = 2.80416; 

chord = S/span; 

mass = 5.9; % 
Ixx = 0.8; % 
Tyy = 0.7; % 
P22 = 1.2; % 


% 


fe) 


& su 
6 wing 
chord 


gross weight 
main moment of inertia around axis Ox (kg*sq.m) 
main moment of inertia around axis Oy (kg*sq.m) 
main moment of inertia around axis Oz (kg*sq.m) 


& Aerodynamic Derivatives (all 





CLO = 0.2115; 
CLa 2.2403; 
CLa_dot = 2.1963; 
Clq = 7.2479; 
CLDe = 0.6868; % 
CDO 0.0195; 
Apolar = 0.0195; 
Al = -0.0049; 
CYb = -1.8449; 
CYDr = 0.2021; & 


si 


Dashes. i 
Tes ug! 
6 li 
o LAE 
ft due 
%& drag 


%& drag 








%& side 
deforce 


rface area of wing (m2) 
span (m) 
(m) 
(kg) 


per radian) 

oefficient at a = 0 = 0.0003; 
urve slope 

ft due to angle of attack rate 
t due to pitch rate 
to elevator 

coefficient at a = 0 

curve slope (A2) 


force due to sideslip 
due to rudder 
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Cma 
Cma_dot 
Cmq 
CmDe 
Cnb 
Cnp 
Cnr 
CnDa 
CnDr 
CLD£ 
CmDf 
% Standa 


sof 


bare 


ISA_lapse = 


ISA_hmax 
ISA_R 
ISA_g 
ISA_rho0O 
ISA_PO 
ISA_TO 


EA ee 


29 8 Ge SC. ee 


-0.1490; 
=3:.1:-97.1.° 
0.4011; 
0.4698; 
0.0103; 
0.0237; 
-3.4836; 
=3..6L14; 
=-13.2651 


ol? 


Atmosphere 
.0065; 
2000; 
287; 
9.815% 
1.225; 
101325; 
289; 


bE ee 2 9.'¢- 


g 2 O-G 8. CO. S 


% dihedral effect =-0.0132 








% roll damping 

% roll due to yaw rate 
roll control power 
ro due to rudder 











% pitch moment at a 0 =>0.0652 

% pitch moment due to angle of attack 

% pitch moment due to angle of attack rate 
pitch moment due to pitch rate 
pitch control power 


z 














%& weathercock stability = 0.075 
% adverse yaw 
% yaw damping 
% aileron adverse yaw 
% yaw control power 
% Lapse rate (degC/m) 
% Altitude limit (m) 
% Gas Constant (degk*m*m/s/s) 
%& Gravity (m/s/s) 
%& Density at sea level (kg/m/m/m) 
% Sea-level Pressure (N/m/m) 
% Sea-level Temperature (degK) 
SEESEESSEEESSEEEEESSEEEEESESS 
SEEEEESEEEESSESEEESSEEEESESESS 
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A. 


APPENDIX B GRAPHICAL UAV MODELS 


P-10B MODEL 


[30x12 | 


12] 25x12 | 30x12 > 


\] 


8.9 












25 30 





10,12 


Figure 55. P-10B Model 
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B. RASCAL MODEL 








4 
x 
1,90 


ia 3.00 3.90 


Figure 56. MS Visio model for Fuselage 


10,30 
3.20 [coma a 
0 





Figure 57. MS Visio model for Wing and Aileron 
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3.30; 14.70 14.7 3.30 
18.0 18.0 


29.40 


(64.9,0,2.40) 


Figure 58. MS Visio model for horizontal stabilizer and Elevator 





Figure 59. MS Visio model for Vertical Stabilizer and Rudder 
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A. P-10B LINAIR INPUT MODEL 


!'!LinAir Pro 


APPENDIX C: LINAIR INPUT MODEL 


! Input file for LinAir Pro for Windows for Rascal 





! Sref Bref Xref 
7620 254 0 

! 

| alpha beta phat 
3 0 0 

! Control surface deflection ( 

! Aileron Elevator Rudder 

! 0 0 0 


Yref 
0 


ghat 
0 
degree) 


Y &&&EEEEEEEEEEEEEEEEEEEEEEEEEEEE EEE EEE EEE 





















































































































































| Wake Pos Reflect CL File 
1 0 CL.dat 

Y && EEE GEE EEE EGE EEEEE EEE EEE EEE EEEEEEEEEEE 

| xLE yLE ZLE xTE 

! Longitudinal 

| Element 1.0000 Fuselag Verticl 

| XrootLE YrootLE ZrootLE XrootTE 
—32.500 0.0000 -6.0000 22.5000 

| XtipLE YtipLk ZtipLE XtipTE 
-32.5000 0.0000 6.0000 22.500 

! Number of Spanwise panels 
2.0000 

| CDO CD1 CD2 Cmac 
0.0035 0.0000 0.0000 0.0000 

| Element 2.0000 Fuselag Verticl 

| XrootLE YrootLE ZrootLE XrootTE 
22.5000 0.0000 -6.0000 52.0000 

| XtipLE YtipLE ZtipLE XtipTE 
22.5000 0.0000 6.0000 52.000 

| Number of Spanwise panels 
3.0000 

| CDO CD1 CD2 Cmac 
0.0035 0.0000 0.0000 0.0000 

| Element 3.0000 Fuselase Vertical 

| XrootLE YrootLE ZrootLE XrootTE 
52.0000 0.0000 -1.2500 94.5000 

! XtipLE YtipLk ZtipLE XtipTE 
52.0000 0.0000 1.2500 94.500 

! Number of Spanwise panels 
4.0000 

! CDO CD1 CD2 Cmac 
0.0035 0.0000 0.0000 0.0000 


a 


Zref Nelem 

0 22 

rhat Mach 

0 0.1 

Out File lem File 





Forces.dat! 














lement.dat 





























yTE ZTE 

1/3 

YrootTE ZrootTE 
0.0000 -—6.0000 
YtipTE ZtipTE 

0.0000 6.0000 

CLmax 

1.4000 

2/3 

YrootTE ZrootTE 
0.0000 =1.2:2500 
YtipTE ZtipTE 

0.0000 1.2500 

CLmax 

1.4000 

3/3 

YrootTE ZrootTE 
0.0000 -1.2500 
YtipTE ZtipTE 

0.0000 1.2500 

CLmax 

1.4000 





Element 
XrootLE 
-—32.5000 
XtipLE 
-—32.5000 
Number 
2.0000 
CDO 
0.0035 














Element 
XrootLE 
22.5000 
XtipLE 
22.5000 
Number 
3.0000 
CDO 
0.0035 














Element 
XrootLE 
52.0000 
XtipLE 
52.0000 
Number 
4.0000 
CDO 
0.0035 














Element 
XrootLE 
-7.5000 
XtipLE 
-7.5000 
Number 
2.0000 
CDO 
0.0100 














Element 
XrootLE 
-7.5000 
XtipLE 
-7.5000 
Number 
4.0000 
CDO 
0.0100 














Elemen 
XroothL 
-7.5000 
XtipLE 





Ei ct 











































































































4.0000 Fuselase 
YrootLE ZrootLE 
—6.0000 0.0000 
YtipLE ZtipLE 
6.0000 0.0000 
of Spanwise 
ep CD2 
0.0000 0.0000 
5.0000 Fuselase 
YrootLE ZrootLE 
—6.0000 0.0000 
YtipLE ZtipLE 
6.0000 0.0000 
of Spanwise 
CD1 CD2 
0.0000 0.0000 
6.0000 Fuselase 
YrootLE ZrootLE 
-1.2500 0.0000 
YtipLE ZtipLE 
1.2500 0.0000 
of Spanwise 
CD1 CD2 
0.0000 0.0000 
7.0000 Left Wing 
YrootLE ZrootLE 
-126.3371 18.6479 
YtipLE ZtipLE 
-54.7316 16.0347 
of Spanwise 
CD1 CD2 
0.0000 0.0418 
8.0000 Left Wing 
YrootLE ZrootLE 
-54.7316 16.0347 
YtipLE ZtipLE 
—6.0000 6.0000 
of Spanwise 
ep CD2 
0.0000 0.0418 
9.0000 Right Wing 
YrootLE ZrootLE 
6.0000 6.0000 
YtipLkE ZtipLE 
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Horizontall/3 

















































































































XrootTE YrootTE ZrootTE 
22.5000 -—6.0000 0.0000 
XtipTE YtipTE ZtipTE 
22.5000 6.0000 0.0000 
panels 
Cmac CLmax 
0.0000 1.4000 
Horizontal2/3 
XrootTE YrootTE ZrootTE 
52.0000 -1.2500 0.0000 
XtipTE YtipTE ZtipTE 
52.0000 1.2500 0.0000 
panels 
Cmac CLmax 
0.0000 1.4000 
Horizontal3/3 
XrootTE YrootTE ZrootTE 
94.5000 -1.2500 0.0000 
XtipTE YtipTE ZtipTE 
94.5000 1.2500 0.0000 
panels 
Cmac CLmax 
0.0000 1.4000 

1/2 
XrootTE YrootTE ZrootTE 
14.0000 -126.3371 18.6479 
XtipTE YtipTE ZtipTE 
14.0000 -54.7316 16.0347 
panels 
Cmac CLmax 
0.0000 1.4000 

2/2 
XrootTE YrootTE ZrootTE 
22.5000 -—54.7316 16.0347 
XtipTE YtipTE ZtipTE 
22.5000 -—6.0000 6.0000 
panels 
Cmac CLmax 
0.0000 1.4000 

1/2 
XrootTE YrootTE ZrootTE 
22.5000 6.0000 6.0000 
XtipTE YtipTE ZtipTE 











-7.5000 
Number 
4.0000 
CDO 
0.0100 





Element 
XrootLE 
-7.5000 
XtipLE 
-7.5000 
Number 
2.0000 
CDO 
0.0100 














Element 
XrootLE 
14.0000 
XtipLE 
14.0000 
Number 
2.0000 
CDO 
0.0100 














Element 
XrootLE 
14.0000 
XtipLE 
14.0000 
Number 
2.0000 
CDO 
0.0100 














Element 
XrootLE 
80.5000 
XtipLE 
80.5000 
Number 
4.0000 
CDO 
0.0090 














Element 
XrootLE 
88.1250 
XtipLE 
88.1250 
Number 
4.0000 
CDO 
0.0090 




















































































































































































































54.7316 16.0347 22.5000 54.7316 16.0347 
of Spanwise panels 

CD1 CD2 Cmac CLmax 

0.0000 0.0418 0.0000 1.4000 

10.0000 Right Wing 2/2 

YrootLE ZrootLE XrootTE YrootTE ZrootTE 
54.7316 16.0347 14.0000 54.7316 16.0347 
YtipLE ZtipLE XtipTE YtipTE ZtipTE 

126.3371 18.6479 14.0000 126.3371 18.6479 
of Spanwise panels 

CD1 CD2 Cmac CLmax 

0.0000 0.0418 0.0000 1.4000 

11.0000 Left Aileron 

YrootLE ZrootLE XrootTE YrootTE ZrootTE 
-126.3371 18.6479 22.5000 -126.3371 18.6479 
YtipLE ZtipLE XtipTE YtipTE ZtipTE 

-—54.7316 16.0347 22.5000 -—54.7316 16.0347 
of Spanwise panels 

CD1 CD2 Cmac CLmax 

0.0000 0.0418 0.0000 1.4000 

12.0000 Right Aileron 

YrootLE ZrootLE XrootTE YrootTE ZrootTE 
54.7316 16.0347 22.5000 54.7316 16.0347 
YtipLE ZtipLE XtipTE YtipTE ZtipTE 

126.3371 18.6479 22.5000 126.3371 18.6479 
of Spanwise panels 

CD1 CD2 Cmac CLmax 

0.0000 0.0418 0.0000 1.4000 

13.0000 Horizontal Stabilizel/2 

YrootLE ZrootLE XrootTE YrootTE ZrootTE 
-—32.5000 1.2500 88.1250 -—32.5000 1.2500 

YtipLE ZtipLE XtipTE YtipTE ZtipTE 

32.5000 1.2500 88.1250 32.5000 1.2500 

of Spanwise panels 

CD1 CD2 Cmac CLmax 

0.0000 0.0706 0.0000 1.4000 

14.0000 Horizontal Stabilize2/2 

YrootLE ZrootLE XrootTE YrootTE ZrootTE 
-—26.1250 1.2500 94.7500 —-26.1250 1.2500 

YtipLE ZtipLE XtipTE YtipTE ZtipTE 

26.1250 1.2500 94.7500 26.1250 1.2500 

of Spanwise panels 

CD1 CD2 Cmac CLmax 

0.0000 0.0706 0.0000 1.4000 
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Element 
XrootLE 
88.1250 
XtipLE 
88.1250 
Number 
2.0000 
CDO 
0.0090 














Element 
XrootLE 
94.7500 
XtipLE 
94.7500 
Number 
2.0000 
CDO 
0.0090 














Element 
XrootLE 
94.7500 
XtipLE 
94.7500 
Number 
2.0000 
CDO 
0.0090 














Element 
XrootLE 
88.1250 
XtipLE 
88.1250 
Number 
2.0000 
CDO 
0.0090 














Element 
XrootLE 
80.5000 
XtipLE 
84.3657 
Number 
4.0000 
CDO 
0.0090 














Element 
XrootLE 
84.3657 


XtipLE 





























































































































































































































15.0000 Left Elevator 1/2 

YrootLE ZrootLE XrootTE YrootTE ZrootTE 
-—32.5000 1.2500 101.7500 -32.5000 1.2500 
YtipLk ZtipLE XtipTE YtipTE ] 
-26.1250 1.2500 101.7500 -26.1250 1.2500 
of Spanwise panels 

ep CD2 Cmac CLmax 

0.0000 0.0706 0.0000 1.4000 

16.0000 Left Elevator 2/2 

YrootLE ZrootLE XrootTE YrootTE ZrootTE 
-26.1250 1.2500 101.7500 -26.1250 1.2500 
YtipLE ZtipLE XtipTE YtipTE ZtipTE 
0.0000 1.2500 101.7500 -4.5000 1.2500 
of Spanwise panels 

CD1 CD2 Cmac CLmax 

0.0000 0.0706 0.0000 1.4000 

17.0000 Right Elevator 1/2 

YrootLE ZrootLE XrootTE YrootTE ZrootTE 
0.0000 1.2500 101.7500 4.5000 1.2500 
YtipLE ZtipLE XtipTE YtipTE ZtipTE 
26.1250 1.2500 101.7500 26.1250 1.2500 
of Spanwise panels 

CD1 CD2 Cmac CLmax 

0.0000 0.0706 0.0000 1.4000 

18.0000 Right Elevator 2/2 

YrootLE ZrootLE XrootTE YrootTE ZrootTE 
26.1250 1.2500 101.7500 26.1250 1.2500 
YtipLE ZtipLE XtipTE YtipTE ZtipTE 
32.5000 1.2500 101.7500 32.5000 1.2500 
of Spanwise panels 

CD1 CD2 Cmac CLmax 

0.0000 0.0706 0.0000 1.4000 

19.0000 Vertical Stabilizer 1/2 

YrootLE ZrootLE XrootTE YrootTE ZrootTE 
0.0000 1.2500 94.7500 0.0000 1.2500 
YtipLE ZtipLE XtipTE ZtipTE ZtipTE 
0.0000 22.1250 94.7500 0.0000 222-4250 
of Spanwise panels 

ep CD2 Cmac CLmax 

0.0000 0.0706 0.0000 1.4000 

20.0000 Vertical Stabilizer 2/2 

YrootLE ZrootLE XrootTE YrootTE ZrootTE 
0.0000 Z2 0025.0 88.8657 0.0000 221250 
YtipLE ZtipLE XtipTE YtipTE ZtipTE 
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85.5000 
! Number 

2.0000 
! CDO 

0.0090 





! Element 
! XrootLE 
94.7500 
! XtipLE 
94.7500 
! Number 
4.0000 
! CDO 
0.0090 











20.13426 

! XrootLE 
88.8657 

! XtipLE 
88.8657 

! Number 
4.0000 

! CDO 
0.0090 


Gp 
G 











end 









































0.0000 28.2500 
of Spanwise 
CD1 CD2 
0.0000 0.0706 
21.0000 Rudder 
YrootLE ZrootLE 
0.0000 1.2500 
YtipLE ZtipLE 
0.0000 22.1250 
of Spanwise 
CD1 CD2 
0.0000 0.0706 
Element 22.0000 
YrootLE ZrootLE 
0.0000 22.1250 
YtipLE ZtipLE 
0.0000 28.2500 
of Spanwise 
CD1 CD2 
0.0000 0.0706 
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88.8657 0.0000 28.2500 
panels 
Cmac CLmax 
0.0000 1.4000 
1/2 

XrootTE YrootTE ZrootTE 
108.7500 0.0000 1.2500 
XtipTE YtipTE ZtipTE 
108.7500 0.0000 22.1250 
panels 
Cmac CLmax 
0.0000 1.4000 

Rudder 108.7500 
XrootTE YrootTE ZrootTE 
108.7500 0.0000 22.1250 
XtipTE YtipTE ZtipTE 
108.7500 0.0000 28.2500 
panels 
Cmac CLmax 
0.0000 1.4000 


2/2 


B. RASCAL LINAIR MODEL 


!'!LinAir Pro 


! Input file for LinAir Pro for Windows for Rascal 























! Sref Bref Xref Yref Zref 
ALD ae TS = AO A 0 0 0 

! 

| alpha beta phat ghat rhat 
3 0 0 0 0 

Y && E&OE EEE EE EE EE EE EE EEE EEE EEE EE KEEEE 

! Wake PosReflect CL File Out File 
tl 0 CL.dat Forces.dat 

Y &&&&EEEEE EEE EE EE EE EE EE EEE EE EEEEEEEEEEE 

| xLE yLE ZLE xTE yTE ZTE 

| Element 1 Fuselage Verticl 1/6 

! XrootLE YrootLE ZrootLE XrootTE 
=23::5 0 0 -20 

! XtipLE YtipLE ZtipLE XtipTE 
-23.5 0 0 -20 

| Number of Spanwise panels 
2 

! CDO CD1 CD2 Cmac 
0.00352 0 0 0 

| Element 2 Fuselage Verticl 2/6 

! XrootLE YrootLE ZrootLE XrootTE 
-20 0 3 —9..97 

! XtipLE YtipLE ZtipLE XtipTE 
-20 0 -3 S957 

| Number of Spanwise panels 
3 

! CDO CD1 CD2 Cmac 
0.00352 0 0 0 

! Element 3 Fuselase Vertical 3/6 

! XrootLE YrootLE ZrootLE XrootTE 
=9..°7 0 4 ah 2 

! XtipLE YtipLE ZtipLE XtipTE 
= O.ey, 0 -4 =132 

| Number of Spanwise panels 
4 

! CDO EDL CD2 Cmac 
0.00352 0 0 0 

! Element 4 Fuselase Vertical 4/6 

! XrootLE YrootLE ZrootLE XrootTE 
Sle2 0 6.7 10.8 

! XtipLE YtipLE ZtipLE XtipTE 
a1 2 0 -4 10.8 

| Number of Spanwise panels 
5 

! CDO CD1 CD2 Cmac 


76B 


Mach 
0.1 


lem File 





7] 











ilement.dat 










































































YrootTE ZrootTE 
0 19 
YtipTE ZtipTE 
0 —129 
CLmax 

1.4 

YrootTE ZrootTE 
0 4 
YtipTE ZtipTk 
0 -4 
CLmax 

1.4 

YrootTE ZroottTE 
0 6.7 
YtipTE ZtipTE 
0 -4 
CLmax 

1.4 

YrootTE ZrootTE 
0 6.7 
YtipTE ZtipTE 
0 -4 
CLmax 


0.00352 


Element 
XrootLE 
10.8 
XtipLE 
10.8 
Number 
3 

CDO 
0.00352 


Element 
XrootLE 
39.9 
XtipLE 
BO 9 
Number 
7 

CDO 
0.00352 


Element 
XrootLE 
=23.49 
XtipLE 
=2335 
Number 
2 

CDO 
0.00352 


Element 
XrootLE 
-20 
XtipLE 
-20 
Number 
3 

CDO 
0.00352 


Element 
XrootLE 
=9..-7 
XtipLE 
-9.7 
Number 
4 

CDO 
0.00352 


Element 
XrootLE 
ee? 


0 0 

5 Fuselase Vertical 5/6 
YrootLE ZrootLE 
0 6.7 
YtipLE ZtipLE 
0 —4 
of Spanwise 
CD1 CD2 
0 0 

6 Fuselase Vertical 6/6 
YrootLE ZrootLE 
0 2.4 
YtipLE ZtipLE 
0 = 2. iil: 
of Spanwise 
CDi CD2 
0 0 

7 Fuselase Horizontal 
YrootLE ZrootLE 
0 0 
YtipLE ZtipLE 
0 0 
of Spanwise 
ED CD2 
0 0 


8 Fuselase 
YrootLE 


CD1 


9 Fuselase 
YrootLE 


CD1 


10 Fuselase 
YrootLE 
4 


Horizontal 


XrootTE 


309 


XtipTE 


39.9 
panel 


Cmac 
0 





LS 


XrootTE 


45.1 
Xtipl 
45.1 
panel 


Cmac 
0 


1/5 


TE 





LS 


XrootTE 


—20 
Xtipl 
—20 
panel 


Cmac 
0 


2/5 





ZrootLE 
0 

ZtipLE 

0 
Spanwise 


CD2 
0 


Horizontal 
ZrootLE 
0 
ZtipLE 
0 
Spanwise 


CD2 
0 


Horizontal 


ZrootLE 
0 


77 


TE 





LS 


XrootTE 


=D] 
Xtipl 
-9.7 
panel 


Cmac 
0 


3/5 


TE 





LS 


XrootTE 


=14.2 
Xtipl 
“1.2 
panel 


Cmac 
0 


4/5 


TE 





LS 


XrootTE 


10.8 





































































































1.4 

YrootTE ZrootTE 
0 2.4 
YtipTE ZtipTk 
0 aed! 
CLmax 

1.4 

YrootTE ZrootTE 
0 2.4 
YtipTE ZtipTE 
0 A cek! 
CLmax 

1.4 

YrootTE ZrootTE 
Tian: 0 
YtipTE ZtipTk 
et eee) 0 

CLmax 

1.4 

YrootTE ZrootTE 
4 0 
YtipTE ZtipTE 
-4 0 

CLmax 

1.4 

YrootTE ZrootTE 
4 0 
YtipTE ZtipTE 
-4 0 

CLmax 

1.4 

YrootTE ZrootTE 
4 0 


XtipLE 
=—1,2 
Number 
5 

CDO 
0.00352 


Element 
XrootLE 
10.8 
XtipLE 
10.8 
Number 
iS! 

CDO 
0.00352 


Element 
XrootLE 
=122 


Element 
XrootLE 
= Ded 
XtipLE 
ee 
Number 

4 

CDO 
0.000033 


Element 
XrootLE 
=1.2 
XtipLE 
-1.2 
Number 

4 

CDO 
0.000033 


Element 
XrootLE 
Ded 
XtipLE 
-3.4 
Number 
4 

CDO 


CD1 


11 Fuselase 
YrootLE 
4 
YtipLE 
-4 
of 


CD1 
0 


12 Right Wing 
YrootLE 


ZtipLE 
0 
Spanwise 


CD2 
0 


Horizontal 
ZrootLE 
0 
ZtipLE 
0 
Spanwise 


CD2 
0 


ZrootLkE 
6.9 
ZtipLE 
T2535 
Spanwise 


CD2 


0 0.020668746 


13 Right Wing 
YrootLE 
-2.8 
YtipLE 
-46 
of 


CD1 
0.0002 


14 Right Wing 
YrootLE 
-21 
YtipLE 
-46 
of 


CD1 
0.0002 


15 Right Wing 
YrootLE 
-46 
YtipLE 
#5932 
of 


CD1 


ZrootLE 
Tegal 
ZtipLE 
8.4 
Spanwise 


CD2 
0.0043 


ZrootLE 
125 
ZtipLE 
8.3 
Spanwise 


CD2 
0.0043 


ZrootLkE 
8.4 
ZtipLE 
8.65 
Spanwise 


CD2 
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XtipTE 





panels 


5/5 


XrootTE 


45.1 


XtipTE 


45.1 


panels 


Cmac 
0 


1/4 


XrootTE 


11 





2/4 


TE 


Is 


XrootTE 


eel ae 
Xtipl 
=? 
panel 





Cmac 
0 


3/4 


TE 


LS 


XrootTE 


8.5 
Xtipl 
4.3 
panel 





Cmac 
0 


4/4 


TE 


LS 


XrootTE 


6.8 
Xtipl 
1.8 
panel 


Cmac 


TE 





LS 





































































































YtipTE ZtipTE 
-4 0 

CLmax 

1.4 

YrootTE ZrootTE 
0.6 0 
YtipTE ZtipTE 
-0.6 0 

CLmax 

1.4 

YrootTE ZrootTE 
0 6.7 
YtipTE ZtipTE 
-21 Presi: 
CLmax 

1.4 

YrootTE ZrootTE 
-2.8 7 
YtipTE ZtipTE 
-46 B33 
CLmax 

1.4 

YrootTE ZrootTE 
-21 7.4 
YtipTE ZtipTE 
-46 8.23 
CLmax 

1.4 

YrootTE ZrootTE 
-46 8.2 
YtipTE ZtipTE 
S042 B50 
CLmax 


0.000033 0.0002 


Element 16 Left Wing 


XrootLE YrootLE 
a2 21 
XtipLE YtipLE 
=4,.2 0 
Number of 

3 

CDO CD1 
0.000033 0.0002 


Element 17 Left Wing 


XrootLE YrootLE 
OZ 46 
XtipLE YtipLE 
252 2.8 
Number of 

4 

CDO CD1 
0.000033 0.0002 


Element 18 Left Wing 


XrootLE YrootLE 
=1.2 46 
XtipLE YtipLE 
-1.2 21 
Number of 

4 

CDO CD1 
0.000033 0.0002 


Element 19 Left Wing 


XrootLE YrootLE 
-3.4 ey 
XtipLE YtipLE 
=92 46 
Number of 

4 

CDO EDs 
0.000033 0.0002 


Element 20 Right Aileron 


XrootLE YrootLE 
8.5 =2 1 
XtipLE YtipLE 
4.3 -46 
Number of 

4 

CDO CD1 
0.000033 0.0002 


Element 21 Left Aileron 


XrootLE YrootLE 
4.3 46 


0.0043 


ZrootLkE 
Leos 
ZtipLE 
6.9 
Spanwise 


CD2 
0.0043 


ZrootLkE 
8.4 
ZtipLE 
Tel 
Spanwise 


CD2 
0.0043 


ZrootLkE 
8.3 
ZtipLE 
Wesyo) 
Spanwise 


CD2 
0.0043 


ZrootLE 
8.65 
ZtipLE 
8.4 
Spanwise 


CD2 
0.0043 


ZrootLE 
7.4 
ZtipLE 
8.23 
Spanwise 


CD2 
0.0043 


ZrootLE 
8.23 
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1/4 


XrootTE 


11 


XtipTE 





panels 


2/4 


XrootTE 


=12 


XtipTE 


-1.2 





panels 


Cmac 
0 


3/4 


XrootTE 


4.3 


XtipTE 


8.5 





panels 


Cmac 
0 


4/4 


XrootTE 


1.8 
Xtipl 
6.8 
panel 


Cmac 
0 


TE 





LS 


XrootTE 


11 
Xtipl 
6.8 
panel 


Cmac 
0 


TE 





LS 


XrootTE 


6.8 





































































































1.4 

YrootTE ZrootTE 
21 hs. 3D 
YtipTE ZtipTE 
0 6.7 
CLmax 

1.4 

YrootTE ZrootTE 
46 843 
YtipTE ZtipTE 
208 7 

CLmax 

1.4 

YrootTE ZrootTE 
46 823 
YtipTE ZtipTk 
21 7.4 
CLmax 

1.4 

YrootTE ZrootTE 
Dose 8.57 
YtipTE ZtipTE 
46 8.2 
CLmax 

1.4 

YrootTE ZrootTE 
-21 rene) 
ZtipTE ZtipTE 
-46 8.2 
CLmax 

1.4 

YrootTE ZrootTE 
46 8.2 


XtipLE 
8.3.5 
Number 

4 

CDO 
0.000033 


Element 22 Horizontal 


XrootLE 
39°.9 
XtipLE 
39:9 
Number 
4 

CDO 
0.009 


Element 23 Horizontal 


XrootLE 
39:29 
XtipLE 
39:9 
Number 
4 

CDO 
0.009 


Element 24 Horizontal 


XrootLE 
39°29 
XtipLE 
39:9 
Number 
4 

CDO 
0.009 


YtipLE 
21 
of 


CD1 
0.0002 


ZtipLE 
7.4 
Spanwise 


CD2 
0.0043 


Stabilizer 





YrootLE 
-14.7 
YtipLE 
14.7 

of 


CD1 
0 


ZrootLE 
2.4 
ZtipLE 
2.4 
Spanwise 


CD2 
0.070637286 


Stabilizer 





YrootLE 
-18 
YtipLE 
= 4) 
of 


CD1 
0 


ZrootLE 
2.4 
Spanwise 


CD2 
0.070637286 


Stabilizer 





YrootLE 
14.7 
YtipLE 
18 

of 


CD1 
0 


Element 25 Left Elevator 


XrootLE 
45.1 
XtipLE 
45.1 
Number 
4 

CDO 
0.009 


YrootLE 
14.7 
YtipLE 
0.5 

of 


CD1 


Element 26 Left Elevator 


XrootLE 
46.1 
XtipLE 
46.1 
Number 
4 

CDO 


YrootLE 
14.7 
YtipLE 
0.5 

of 


CD1 


ZrootLE 
2.4 
ZtipLE 
2.4 
Spanwise 


CD2 
0.070637286 


ZrootLE 
2.4 
ZtipLE 
2.4 
Spanwise 


CD2 
0.070637286 


ZrootLE 
2.4 
ZtipLE 
2.4 
Spanwise 


CD2 


80 


XtipTE 





panels 


1/3 


XrootTE 


45.1 


XtipTE 


45.1 


panels 


Cmac 
0 


2/3 


XrootTE 


44.1 
Xtipl 
46.1 
panel 





Cmac 
0 


3/3 


TE 


LS 


XrootTE 


46.1 
Xtipl 
44.1 
panel 





Cmac 
0 


1/2 


TE 


LS 


XrootTE 


46.1 
Xtipl 
46.1 
panel 





Cmac 
0 


2/2 


TE 


LS 


XrootTE 


48.3 
Xtipl 
48.3 
panel 


Cmac 


TE 





LS 





































































































YtipTE ZtipTE 
21 Paresrs) 
CLmax 

1.4 

YrootTE ZrootTE 
-14.7 2.4 
YtipTE ZtipTE 
14.7 24 
CLmax 

1.4 

YrootTE ZrootTE 
-18 2.4 
YtipTE ZtipTE 
-14.7 24 
CLmax 

1.4 

YrootTE ZrootTE 
14.7 2.4 
YtipTE ZtipTE 
18 24 
CLmax 

1.4 

YrootTE ZrootTE 
14.7 2.4 
YtipTE ZtipTE 
Oi5 2.4 
CLmax 

1.4 

YrootTE ZroottTE 
8.2 2.4 
YtipTE ZtipTE 
3.7 2.4 
CLmax 


0.009 0) 0.070637286 


Element 27 Right Elevator 


XrootLkE YrootLE ZrootLE 
45.1 =0..5 2.4 

XtipLE YtipLE ZtipLE 

45.1 -14.7 2.4 

Number of Spanwise 

4 

CDO CD1 CD2 

0.009 0 0.070637286 


Element 28 Right Elevator 


XrootLkE YrootLE ZrootLE 
46.1 =0.25 2.4 

XtipLE YtipLE ZtipLE 

46.1 -14.7 2.4 

Number of Spanwise 

4 

CDO CD1 CD2 

0.009 0 0.070637286 


Element 29 Vertical Stabilizer 


XrootLkE YrootLE ZrootLE 
36.1 0 2.4 

XtipLE YtipLE ZtipLE 

41.6 0 11.4 

Number of Spanwise 

4 

CDO CD1 CD2 

0.009 0 0.070637286 


Element 30 Vertical Stabilizer 


XrootLkE YrootLE ZrootLE 
41.6 0 11.4 

XtipLE YtipLE ZtipLE 

43.6 0 134:6 

Number of Spanwise 

4 

CDO CD1 CD2 

0.009 0 0.070637286 


Element 31 Rudder 


XrootLkE YrootLE ZrootLE 
45.1 0 2.4 

XtipLE YtipLE ZtipLE 

45.1 0 11.4 

Number of Spanwise 

4 

CDO CD1 CD2 

0.009 0 0.070637286 


Element 32 Rudder 
XrootLE YrootLE ZrootLE 
45.1 0 -0.6 
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Uf? 
XrootTE 
46.1 
XtipTE 
46.1 
panels 





Cmac 
0 


2/2 
XrootTE 
438.3 
XtipTE 
48.3 
panels 





Cmac 
0 


1/2 
XrootTE 
45.1 
XtipTE 
45.1 
panels 





Cmac 
0 


2/2 
XrootTE 
47.3 
XtipTE 
45.3 
panels 





Cmac 
0 


1/2 
XrootTE 
49.6 
XtipTE 
47.3 
panels 





Cmac 
0 


2/2 
XrootTE 
46.6 





































































































1.4 

YrootTE ZrootTE 
=015 2.4 
YtipTE ZtipTk 
-14.7 2.4 
CLmax 

1.4 

YrootTE ZrootTE 
3.5 1 2.4 
YtipTE ZtipTE 
-8.2 2.4 
CLmax 

1.4 

YrootTE ZrootTE 
0 2.4 
ZtipTeE ZtipTeE 
0 11.4 
CLmax 

1.4 

YrootTE ZrootTE 
0 11.4 
YtipTE ZtipTE 
0 1376 
CLmax 

1.4 

YrootTE ZrootTE 
0 2.4 
ZtipTE ZtipTE 
0 11.4 
CLmax 

1.4 

YrootTE ZrootTE 
0 -0.6 


end 


XtipLE 
45.1 
Number 
4 

CDO 
0.009 


ZtipLE 
2.4 
Spanwise 


CD2 
0.070637286 
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XtipTE 
49.6 
panels 





Cmac 
0 











YtipTE ZtipTE 
0 264 
CLmax 

1.4 


APPENDIX D: SIMULINK TRIM COMMAND 


These commend takes input from both the P10B and Rascal inout file to compute 


the trim condition for each UAV. 





6% Determine steady-state trim parameters 

cle, V=norm(Vb_0); df=0; 

ro=1.225* (1-0.0065*Pos_0(3) /288.15) *4.25589; 

atrim = (CmDe*CLO*S*V*%2*rot+tCmDe*CLDE*df*S*V*%2*ro-2*CmDe*mass*ISA_g- 
CLDe*S*V*2*ro*Cm0-CLDe*S*V*2*ro*CmDf*df) /S/ro/V*2/ (-CLa*CmDe+Cma*CLDe) ; 


det rim= — (Cma*CLO*S*V*%2*ro+Cma*CLDf*df*S*V*2*ro-2*Cma*mass*ISA_g- 
CLa*S*V*2*ro*Cm0-CLa*S*V*%2*ro*CmD£*df) /S/ro/V*2/ (-CLa*CmDe+Cma*CLDe) ; 





atrim=atrim*180/pi 
det rim=detrim*180/pi 
CLtrim=mass*ISA_g/ (S*ro*V*2/2) ; 


Thrusttrim=(CDO+A1*CLtrim+tApolar*CLtrim’2) *S*ro*V%2/2 
dthrottle=Thrusttrim/244.5 





83 


THIS PAGE INTENTIONALLY LEFT BLANK 


84 


APPENDIX E: AVL INPUT FORMAT 


A sample file for an RC glider is shown below. Comment lines begin with a "#". 


myn 


Everything after and including a "!" is ignored. Blank lines are ignored. 

# SuperGee 

# 

# Dimensional unit and parameter data. 

# Mass & Inertia breakdown. 

# Names and scalings for units to be used for trim and eigenmode calculations. 
# The Lunit and Munit values scale the mass, xyz, and inertia table data below. 
# Lunit value will also scale all lengths and areas in the AVL input file. 

Lunit = 0.0254 m 

Munit = 0.001 kg 


Tunit = 1.0s 


# Gravity and density to be used as default values in trim setup (saves runtime 
typing). 


# Must be in the unit names given above (i.e. m,kg,s). 


g =9.81 
rho = 1.225 
ge aire eh eee: Beier eee 


# Mass & Inertia breakdown. 


# x y Z 1s location of item's own CG. 
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# Ixx... are item's inertias about item's own CG. 
# 
# X,y,z system here must be exactly the same one used in the .avl input file 


# (same orientation, same origin location, same length units) 


# mass x y Z_ Ixx Ilyy Izz [ Ixy Ixz lyz] 

sain Commas Cues Nama Cea ous Mem I, he, 

+0 0 0 0 0 0 0. 0 0 0. 

58.0 3.34 12.0 1.05 4400 180 4580 ! right wing 
58.0 3.34-12.0 1.05 4400 180 4580 ! left wing 
16.0 -5.2 0.0 0.0 O 80 80 ! fuselage pod 
18.0 13.25 0.0 0.0 0 700 700 ! boom+rods 
22.0 -7.4 00 0.0 O O 0 ! battery 
2.0 -25 0.0 0.0 O O 0 ! jack 
9.0 -3.8 0.0 0.0 O O 0 !'RX 
9.0 -5.1 00 00 O O 0 ! rud servo 
6.0 -5.9 0.0 0.0 O O 0 ! ele servo 
90 26 10 00 0 O 0 ! R wing servo 
9.0 26-10 00 O O 0 ! L wing servo 
20 1000 05 0 0O 0 ! wing connector 
10 30 00 00 0 0O 0 ! wing pins 
6.0 29.0 0.0 10 70 2 72 ! stab 
6.0 33.0 0.0 2.0 35 39 4 ! rudder 


0.0 -8.3 00 00 0 0O 0 ! nose wt. 
86 


Units 


The first three lines give the magnitudes and names of the units to be used for run 


case setup and possibly for eigenmode calculations. 
o Lunit = 0.0254 m 
o Munit=0.001 kg 
o Tunit=1.0s 


In this example, standard SI units (m,kg,s) are chosen. But the data in xxx.avl and 
XXX.mass is given in units of Lunit = 1 inch, which is therefore declared here to be equal 
to "0.0254 m". If the data was given in centimeters, the statement would read, “Lunit = 


0.01 m”’, and if it was given directly in meters, it would read, “Lunit = 1.0 m” 


Similarly, Munit used here in this file is the gram, but since the kilogram (kg) is to 
be used for run case calculations, the Munit declaration is, “Munit = 0.001 kg”. If the 
masses here were given in ounces, the declaration would be, “Munit = 0.02835 kg”. The 
third line gives the time unit name and magnitude. If any of the three unit lines is absent, 
that unit's magnitude will be set to 1.0, and the unit name will simply remain as "Lunit", 


"Munit", or "Tunit" 
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APPENDIX F: AVL INPUT FILE 


A. P-10B MODEL 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KKK KKKKK 


AVL dataset for Untitled Aircraft model 
Generated by AVL Model Editor on 28 Oct 2008 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 


Untitled Aircraft 

















Mach 

0.0920 

TYsym IZsym Zsym 

0 0 0.0000 
Sref Cref Bref 
@Auto-generate 
1815.000015.0000 121.0000 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KKK KKKKK 





AVL Axes: 
+X downstream 
EY: out right wing 
+Z up 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KKK KKKKK 





Xref Yref Zref 
0.0000 0.0000 2.5775 
CDp 

0.0000 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KKK KKKKK 


Surfaces 
KKKKKKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK K KK KK KKK KKK 























Wing 

SURFACE 
Wing 

Nchord Cspace Nspan Sspace 
12 2.0000 64 0.0000 
SCALE 
#SX sy sZ 

1.0000 1.0000 1.0000 
TRANSLATE 

#dX dy dz 

-3.7500 0.0000 3.0000 
ANGLE 
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# Wing section 1 
SECTION 

#XLle Yle Zle Chord Angle 
0.0000 -60.5000 0.0000 15.0000 3.5000 





AFILE 
#Airfoil definition 
clarky.txt 








CLAF 
#CLaf = CLalpha / (2 * pi) 
1.0000 


CONTROL 
#label gain Xhinge xhvec Yhvec Zhvec 
LeftAileron -1.0000 0.7167 0.0000 0.0000 0.0000 


SgnDup 








# Wing section 2 











SECTION 
#Xle Yle Zle Chord Angle 
0.0000 -48.0000 0.0000 15.0000 3.5000 


AFILE 
#Airfoil definition 
clarky.txt 








CLAF 
#CLaf = CLalpha / (2 * pi) 
1.0000 


CONTROL 
#label gain Xhinge xhvec Yhvec Zhvec 
LeftAileron -1.0000 0.7167 0.0000 0.0000 0.0000 


SgnDup 
1 

















# Wing section 3 
SECTION 
#Xle Yle Zle Chord Angle 


0.0000 0.0000 0.0000 15.0000 3.5000 














AFILE 
#Airfoil definition 
clarky.txt 
CLAF 
#CLaf = CLalpha / (2 * pi) 
1.0000 
# Wing section 4 
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SECTION 
#Xle Yle Zle Chord Angle 
0.0000 48.0000 0.0000 15.0000 3.5000 
































AFILE 
#Airfoil definition 
clarky.txt 
CLAF 
#CLaf = CLalpha / (2 * pi) 
1.0000 
CONTROL 
#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
RightAileron -1.0000 0.7167 0.0000 0.0000 0.0000 1 
+ Wing section 5 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 60.5000 0.0000 15.0000 3.5000 
AF ILE 
#Airfoil definition 
clarky.txt 
CLAF 
#CLaf = CLalpha / (2 * pi) 
1.0000 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
RightAileron -1.0000 0.7167 0.0000 0.0000 0.0000 1 




















# Horizontal 
#@Yduplicate 3 0.00000 Horizontal 
SURFACE 
RHorizontal 
#Nchord Cspace Nspan Sspace 

12 2.0000 32 0.0000 
SCALE 
#SX sy sZ 


1.0000 1.0000 1.0000 


TRANSLATI 
#dxX dy dZ 
40.2500 0.0000 0.6250 





| 


ANGLE 
#Ainc 
0.0000 
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#Lsurf 








































































































1 
# Horizontal section 1 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 16.2500 0.0000 10.5000 0.0000 
NACA 
#Airfoil definition 
0008 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Basename Elevator 
RElevator -1.0000 0.3512 14.0000 0.0000 0.0000 dl! 
# Horizontal section 2 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 13.0625 0.0000 10.5000 0.0000 
NACA 
#Airfoil definition 
0008 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Basename Elevator 
RElevator -1.0000 0.6667 0.0000 0.0000 0.0000 -1 
# Horizontal section 3 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 0.0000 0.0000 10.5000 0.0000 
NACA 
#Airfoil definition 
0008 
CONTROL 
#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
#@Basename New 
LNew -1.0000 0.6667 0.6667 0.0000 0.0000 pe 
CONTROL 
#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
#@Ignore 
RNew -1.0000 0.6667 0.6667 0.0000 0.0000 = 
# Horizontal (mirror) 
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#@Ignore 
SURFACE 
LHorizonta 
#Nchord 

12 





SCALI 
#SX 
1.0000 





Gl 


TRANSLATE 
#dxX 
40.2500 





ANGLE 
#Ainc 
0.0000 











INDEX 
#Lsurft 
1 


als 
Cspace Nspan 
2.0000 32 


sy sZ 
1.0000 1.0000 


dy dZ 
0.0000 0.6250 








# 

SECTION 

#Xle 
0.0000 





NACA 
#Airfoil d 


Horizontal 


Sspace 
0.0000 


section 4 











Yle 
0.0000 


Zle 
0.0000 


efinition 





0008 


CONTROL 
#label 
#@Basename 

LNew 





CONTROL 

#label 

#@Ignore 
RNew 








Xhinge 


0.6667 


Xhinge 


0.6667 





# 





SECTION 
#Xle 
0.0000 





NACA 
#Airfoil 


Horizontal 


Chord 
10.5000 


Xhvec 


0.6667 


Xhvec 


0.6667 


section 5 











Yle Zle 
-13.0625 0.0000 


definition 





0008 


CONTROL 
#label 
#@Basename 


gain 
Elevator 


Xhinge 





Chord 
10.5000 


Xhvec 


Angle 

0.0000 

Yhvec Zhvec SgnDup 
0.0000 0.0000 = 
Yhvec Zhvec SgnDup 
0.0000 0.0000 cal 
Angle 

0.0000 

Yhvec Zhvec SgnDup 


LElevator -1.0000 0.6667 0.0000 





0.0000 


0.0000 








# Horizontal section 6 











SECTION 
#Xle Yle Zle 
0.0000 -16.2500 0.0000 





Chord 
10.5000 


Angle 
0.0000 


NACA 
#Airfoil 
0008 





definition 


CONTROL 

#label gain 
#@Basename Elevator 
LElevator —-1.0000 


Xhinge xXhvec Yhvec 











0.3512 14.0000 





0.0000 


Zhvec 


0.0000 


SgnDup 








# Vertical 








#@Yduplicate 3 0.00000 Vertical 

SURFACE 

RVertical 

#Nchord Cspace Nspan 
12 2.0000 32 





e5 


Sspace 
0.0000 


SCALI 
#SX sy sZ 
1.0000 1.0000 1.0000 





Gl 


TRANSLATE 
#dx dy dZ 
40.2500 0.0000 0.0000 





ANGLE 
#Ainc 
0.0000 











INDEX 
#Lsurft 
2 














# Vertical section 1 
SECTION 

#Xle Yle Zle 
0.0000 0.5000 0.6250 





Chord 
14.1250 


Angle 
0.0000 


NACA 
#Airfoil 
0008 


definition 





CONTROL 

#label gain 
#@Basename New 
RNew 0.0000 


Xhinge Xhvec Yhvec 


0.5044 0.0000 0.0000 
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Zhvec 


0.0000 


SgnDup 














# Vertical section 2 
SECTION 

#Xle Yle Zle Chord Angle 
2.1829 0.0000 11.0625 12.1921 0.0000 






























































NACA 
#Airfoil definition 
0008 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Basename New 
LNew 0.0000 0.4932 0.0000 0.0000 0.0000 =a 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Ignore 
RNew 0.0000 0.4932 0.0000 0.0000 0.0000 =a! 
# Vertical section 3 
SECTION 
#X1le Yle Zle Chord Angle 
2.7500 0.0000 14.1250 11.6750 0.0000 
NACA 
#Airfoil definition 
0008 
CONTROL 
#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
#@Basename New 
LNew 0.0000 0.1935 0.0000 0.0000 0.0000 = 
CONTROL 
#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
#@Ignore 
RNew 0.0000 0.1935 0.0000 0.0000 0.0000 =A 
# Vertical (mirror) 
#@Ignore 
SURFACE 
LVertical 
#Nchord Cspace Nspan Sspace 
12 2.0000 B32 0.0000 
SCALE 
#sX sY sZ 


1.0000 1.0000 1.0000 


TRANSLAT! 
#dX ay dz 





Gl 


a5 
















































































40.2500 0.0000 0.0000 
ANGLE 
#Ainc 
0.0000 
INDEX 
#Lsurf 
2 
# Vertical section 4 
SECTION 
#Xle Yle Zle Chord Angle 
2.7500 0.0000 14.1250 11.6750 0.0000 
NACA 
#Airfoil definition 
0008 
CONTROL 
#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
#@Basename New 
LNew 0.0000 0.1935 0.0000 0.0000 0.0000 Sal 
CONTROL 
#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
#@Ignore 
RNew 0.0000 0.1935 0.0000 0.0000 0.0000 = 
# Vertical section 5 
SECTION 
#Xle Yle Zle Chord Angle 
2.1829 0.0000 11.0625 12.1921 0.0000 
NACA 
#Airfoil definition 
0008 
CONTROL 
#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
#@Basename New 
LNew 0.0000 0.4932 0.0000 0.0000 0.0000 al, 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Ignore 
RNew 0.0000 0.4932 0.0000 0.0000 0.0000 -1 
# Vertical section 6 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 -0.5000 0.6250 14.1250 0.0000 
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NACA 





#Airfoil definition 

0008 
CONTROL 
#label gain Xhinge 
#@Basename New 

LNew 0.0000 0.5044 


Bodies 


Xhvec 


0.0000 


Yhvec Zhvec SgnDup 


0.0000 0.0000 =. 


KKEKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 





P10B 




















BODY 
P10B 
Nbody Bspace 
15 2.0000 
SCALE 
#SX sy sZ 
1.0000 1.0000 1.0000 
TRANSLATE 
#dX dy dz 
-22.2500 0.0000 0.0000 
BFILE 
#Body file 
P10B.dat 


oF 


B. RASCAL MODEL 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 


AVL dataset for Rascal model 


Generated by AVL Model Editor on 2 Nov 2008 
KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KK KK 











Rascal 

Mach 

0.0588 

TYsym IZsym Zsym 

0 0 0.0000 
Sref Cref Bref 








@Auto-generate 
1509.440013.6725 110.4000 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KK KKK KK 





AVL Axes: 
+X downstream 
EY, out right wing 
+Z up 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKAKKKKKKKKKKKKKKKKKKK KKK KKKKK 





Xref Yref Zref 
0.0000 0.0000 0.0000 
CDp 

0.0000 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK KK 


Surfaces 
KKKKKKKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK 








HorizontalStabiliser 
@Yduplicate 5 0.00000 HorizontalStabiliser 
































SURFACE 
RHorizontalStabiliser 

Nchord Cspace Nspan Sspace 
12 2.0000 32 0.0000 
SCALE 
#SX sy sZ 


1.0000 1.0000 1.0000 


TRANSLAT] 
#dX dy dz 
41.4000 0.0000 1.5000 





GI 


ANGLE 
#Ainc 
0.0000 
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INDEX 


























































































































#Lsurf 

1 
# HorizontalStabiliser section 1 
SECTION 
#Xle Yle Zle Chord Angle 

0.0000 0.0000 0.0000 6.2000 0.0000 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Basename Elevator 

LElevator 0.0000 0.8387 0.0000 0.0000 0.0000 dl! 
CONTROL 

#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Ignore 

RElevator 0.0000 0.8387 0.0000 0.0000 0.0000 1 

# HorizontalStabiliser section 2 

SECTION 

#Xle Yle Zle Chord Angle 

0.0000 3.7000 0.0000 8.4000 0.0000 
CONTROL 

#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
#@Basename Elevator 

RElevator 0.0000 0.6190 0.0000 0.0000 0.0000 1 

# HorizontalStabiliser section 3 

SECTION 

#Xle Yle Zle Chord Angle 

0.0000 8.2000 0.0000 8.4000 0.0000 
CONTROL 

#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Basename Elevator 

RElevator 0.0000 0.6190 0.0000 0.0000 0.0000 1 

# HorizontalStabiliser section 4 

SECTION 

#Xle Yle Zle Chord Angle 

0.0000 14.7000 0.0000 6.2000 0.0000 

CONTROL 

#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Basename Elevator 

RElevator 0.0000 0.8387 0.0000 0.0000 0.0000 1 

















HorizontalStabiliser section 5 
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SECTION 
#Xle Yle Zle Chord Angle 
0.0000 18.0000 0.0000 4.2000 0.0000 





























# HorizontalStabiliser (mirror) 
#@Ignore 
SURFACE 
LHorizontalStabiliser 
#Nchord Cspace Nspan Sspace 
12 2.0000 32 0.0000 
SCALE 
#SX sy sZ 


1.0000 1.0000 1.0000 





1 | 


TRANSLATI 
#dxX dy dz 
41.4000 0.0000 1.5000 
















































































ANGLE 
#Ainc 
0.0000 
INDEX 
#Lsurf 
1 
# HorizontalStabiliser section 6 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 -18.0000 0.0000 4.2000 0.0000 
# HorizontalStabiliser section 7 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 -14.7000 0.0000 6.2000 0.0000 
CONTROL 
#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
#@Basename Elevator 
LElevator 0.0000 0.8387 0.0000 0.0000 0.0000 1 
# HorizontalStabiliser section 8 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 -8.2000 0.0000 8.4000 0.0000 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Basename Elevator 
LElevator 0.0000 0.6190 0.0000 0.0000 0.0000 1 
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HorizontalStabiliser section 9 





# 


































































































SECTION 
#Xle Yle Zle Chord Angle 
0.0000 -3.7000 0.0000 8.4000 0.0000 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Basename Elevator 
LElevator 0.0000 0.6190 0.0000 0.0000 0.0000 1 
# HorizontalStabiliser section 10 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 0.0000 0.0000 6.2000 0.0000 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Basename Elevator 
LElevator 0.0000 0.8387 0.0000 0.0000 0.0000 1 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
#@Ignore 
RElevator 0.0000 0.8387 0.0000 0.0000 0.0000 1 
# VerticalStabiliser 
SURFACE 
VerticalStabiliser 
#Nchord Cspace Nspan Sspace 
12 2.0000 64 0.0000 
SCALE 
#SX sy sZ 
1.0000 1.0000 1.0000 
TRANSLATE 
#dX dy daz 
46.6000 0.0000 1.5000 
ANGLE 
#Ainc 
0.0000 
# VerticalStabiliser section 1 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 0.0000 —-3.0000 1.5000 0.0000 
CONTROL 
#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
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Rudder 0.0000 0.0001 0.0000 0.0000 0.0000 1 

# VerticalStabiliser section 2 

SECTION 

#Xle Yle Zle Chord Angle 

-—9.0000 0.0000 0.0000 13.5000 0.0000 

CONTROL 

#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
Rudder 0.0000 0.6667 0.0000 0.0000 0.0000 1 

# VerticalStabiliser section 3 

SECTION 

#Xle Yle Zle Chord Angle 

—-3.5000 0.0000 9.0000 5.7000 0.0000 

CONTROL 

#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
Rudder 0.0000 0.6140 0.0000 0.0000 0.0000 1 

















# VerticalStabiliser section 4 
SECTION 
#Xle Yle Zle Chord Angle 


-1.5000 0.0000 11.2000 1.7000 0.0000 














# Wing 

SURFACE 

Wing 

#Nchord Cspace Nspan Sspace 
12 2.0000 64 0.0000 

SCALE 

#SX sy sZ 


1.0000 1.0000 1.0000 


TRANSLATE 
#dxX dy dz 
0.0000 0.0000 5.3500 





ANGLE 
#Ainc 
0.0000 











# Wing section 1 








SECTION 
#Xle Yle Zle Chord Angle 
1.8000 -55.2000 1.7000 5.2000 1.0000 





NACA 
#Airfoil definition 
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2312 










































































# Wing section 2 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 -46.0000 1.4000 11.0000 1.0000 
NACA 
#Airfoil definition 
2312 
CONTROL 
#label gain Xhinge Xhvec Yhvec Zhvec SgnDup 
LeftAileron 0.0000 0.7727 0.0000 0.0000 0.0000 1 
+ Wing section 3 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 -21.0000 0.6500 16.2000 1.0000 
NACA 
#Airfoil definition 
2312 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
LeftAileron 0.0000 0.8457 0.0000 0.0000 0.0000 1 
# Wing section 4 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 0.0000 0.0000 16.2000 1.0000 
NACA 
#Airfoil definition 
2312 
# Wing section 5 
SECTION 
#Xle Yle Zle Chord Angle 
0.0000 21.0000 0.6500 16.2000 1.0000 
NACA 
#Airfoil definition 
2312 
CONTROL 
#label gain Xhinge xXhvec Yhvec Zhvec SgnDup 
RightAileron 0.0000 0.8457 0.0000 0.0000 0.0000 ui 
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# 








SECTION 
#Xle Yle 
0.0000 





NACA 


#Airfoil definition 


46.0000 





2312 


CONTROL 
#label 


gain 
RightAileron 0.0000 


Wing section 6 





Zle Chord Angle 

1.4000 11.0000 1.0000 
Xhinge Xhvec Yhvec Zhvec 
OT 727 0.0000 0.0000 0.0000 





SgnDup 














# Wing section 7 
SECTION 

#Xle Yle Zle Chord Angle 
1.8000 55.2000 1.7000 5.2000 1.0000 

NACA 


#Airfoil definition 





2312 


Bodies 





KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 


KKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKAKKKKKKKKKKKKKKKKKKKKKK KKK KK 




















BODY 
Rascal 
Nbody Bspace 
40 1.5000 
SCALE 
#SX sy 
1.0000 1.0000 
TRANSLATE 
#dX dy 


-20.0000 0.0000 


BFILE 
#Body file 
Rascal.dat 








Rascal 


sZ 
1.0000 


dz 
0.0000 
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